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INTRODUCTION 


Due  to  the  serious  interest  in  radiated  electromagnetic  effects  which  exists  within 
the  Department  of  Defense  and  the  electromagnetic  compatibility  (EMC)  community  the 
subject  of  electromagnetic  shielding  has  risen  to  a  position  of  high  importance.  Military 
systems  and  support  equipment  are  generally  required  to  operate  either  during  or  after 
exposure  to  a  variety  of  intense  electromagnetic  environments,  such  as  electromagnetic 
pulse  (EMP),  electromagnetic  interference  (EMI),  jammers,  and  directed  energy 
weapons,  that  is,  high  power  microwaves  (HPM).  The  problem  of  shielding  by  circular 
conducting  cylinders  is  especially  important  for  the  analysis  of  military  systems  due  to 
the  congruity,  or  similarity  of  this  shape  to  the  shapes  of  airplanes,  missiles,  and  artillery 
projectiles.  Rough  engineering  estimates  of  the  shielding  inherent  in  these  types  of 
military  systems  may  be  made  by  approximating  these  structures  by  ideal  cylindrical 
conducting  shells. 

The  propagation  of  an  electromagnetic  planewave  in  the  presence  of  a  planar 
conducting  sheet,  an  infinitely  long  cylindrical  conducting  surface,  and  a  spherical 
conducting  surface  are  among  the  few  problems  which  admit  an  exact  analytical  solu¬ 
tion  and  have  therefore  been  investigated  in  some  detail.  These  solutions  have  been 
used  to  make  calculations  of  the  electromagnetic  attenuation,  or  shielding  effectiveness, 
of  these  structures.  More  complicated  structures  have  also  been  investigated  by  vari¬ 
ous  methods.  The  shielding  effectiveness  is  the  primary  quantity  of  interest  in  the 
subject  of  electromagnetic  shielding.  It  is  a  quantitative  measure  of  the  attenuation  of 
incident  electromagnetic  energy  due  to  the  presence  of  a  metallic  surface  or  shield.  It 
may  be  roughly  defined  as  the  intensity  of  the  electromagnetic  field  at  a  point  in  space 
within  the  shielded  volume  to  the  intensity  of  the  incident  planewave. 

The  purpose  of  this  report  is  to  calculate  the  shielding  effectiveness  versus  fre¬ 
quency  for  circular  cylindrical  shields  and  for  spherical  shields.  An  important  result  of 
this  paper  is  the  derivation  of  an  approximate  formula  which  expresses  the  shielding 
effectiveness  as  a  function  of  the  frequency,  the  conductivity,  and  the  thickness  of  the 
shield.  The  results  show  that  this  formula  is  a  good  approximation  for  frequencies  in  the 
resonance  regime.  This  approximate  formula  permits  quick  shielding  calculations  to  be 
made  without  the  need  for  computer  programs  or  complicated  mathematics. 

Due  to  its  geometrical  and  mathematical  simplicity,  the  problem  of  scattering  of  a 
planewave  by  a  solid  circular  cylinder  at  perpendicular  incidence  was  solved  as  early  as 
1881  by  Lord  Rayleigh,  (ref  1),  and  independently  by  von  Ignatowski  in  1905  (ref  2). 
Following  these  original  works  many  German  authors  investigated  this  problem,  includ¬ 
ing  von  Ignatowski  in  1907,  (ref  3),  Seitz  in  1905  and  1906,  (refs  4-6),  von  Schaefer  in 
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1907  (refs  7,8)  and  Debye  in  1908  (ref  9).  The  scattering  of  a  normally  incident 
planewave  by  a  solid  cylinder  with  a  concentric  cylindrical  sheath  was  first  solved  by 
Adey  in  1956  (ref  10),  Tang  in  1957  (ref  11)  and  latter  investigated  by  Kerker  and 
Matijevic  in  1961  (ref  12).  The  problem  of  scattering  by  a  solid  cylinder  due  to  a 
planewave  at  oblique  incidence  was  not  solved  for  many  years  following  the  initial  work 
of  Lord  Rayleigh  and  von  Ignatowski.  The  full  solution  of  this  problem  was  given  by 
Wait  in  1955  (ref  13)  although  the  method  and  form  of  the  solution  was  given  in  the 
classic  textbook  on  electromagnetic  theory  by  Stratton  in  1941  (ref  14).  Exact  solutions 
for  scattering  by  cylinders  due  to  other  types  of  sources,  such  as  line  sources  and 
dipole  sources,  may  be  found  in  the  encyclopedic  reference  book  on  electromagnetic 
scattering  by  Bowman,  Senior,  and  Uslenghi  (ref  1 5). 

The  first  analysis  of  the  shielding  properties  of  hollow  conducting  cylinders  was 
made  by  Wu  and  Tsai  in  1974  (ref  16).  These  authors  solved  the  scattering  problem  for 
a  cylindrical  conducting  shell  of  inner  radius  r=a  and  outer  radius  r=b  due  to  a  normally 
incident  planewave,  and  then  calculated  the  shielding  effectiveness  as  a  function  of  ka 
for  a  cylinder  of  radius  a=3.0  mm  and  thickness  b-a=0.1  mm.  The  calculation  covered 
the  frequency  range  from  ka=1.1  to  ka=4.0,  and,  therefore,  was  limited  to  a  narrow  band 
in  the  vicinity  of  the  first  resonant  frequency.  In  a  second  paper,  Wu  and  Tsai  analyzed 
the  shielding  effectiveness  for  a  planewave  at  oblique  incidence  (ref  17),  and  in  a  third 
paper  they  investigated  the  transient  response  of  a  cylindrical  shield  to  simple  EMP 
wave  form  (ref  18).  These  are  the  only  results  on  the  shielding  problem  for  circular 
conducting  cylinders  which  are  known  to  this  author.  In  this  report  some  of  this  earlier 
work  shall  be  extended.  The  shielding  effectiveness  shall  be  calculated  over  a  wide 
range  of  frequencies  which  covers  several  orders  of  magnitude,  in  reference  16  no 
details  of  the  numerical  calculations  were  given  other  than  a  reference  to  a  well  known 
handbook  of  functions  (ref  19).  In  this  report,  full  details  of  the  numerical  calculations 
shall  be  included  for  completeness.  This  analysis  leads  to  a  derivation  of  the  resonance 
condition  and  a  derivation  of  an  approximate  formula  for  the  shielding  effectiveness  in 
the  resonance  regime.  No  results  for  the  shielding  effectiveness  are  presented  for  the 
case  of  an  obliquely  incident  planewave  due  to  time  constraints  on  this  work. 

In  the  second  part  of  this  study  the  shielding  properties  of  a  spherical  shield  are 
investigated.  This  problem  is  closely  related  to  the  problem  for  the  cylinder  and  the 
results  are  very  similar.  The  complete  theory  for  scattering  of  an  electromagnetic 
planewave  by  a  sphere  is  usually  attributed  to  Mie  (ref  20).  However,  as  pointed  out  by 
Kerker  (ref  21)  the  theory  had  been  worked  out  earlier  by  several  authors,  and  even  as 
early  as  1863  by  Clebsch.  Some  of  the  earlier  work  included  the  theory  of  scattering  by 
a  perfectly  conducting  sphere  by  Lorenz  in  1 890,  the  theory  of  scattering  by  a  perfectly 
conducting  sphere  by  J.J.  Thomson  in  1893,  and  an  extension  of  Thomson’s  theory  to 
dielectric  and  imperfectly  conducting  spheres  by  Love  in  1899.  Following  the  work  of 
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Mie  in  1908,  further  contributions  to  the  theory  were  made  by  Debye  in  1909  (ref  22). 
The  most  concise  formulations  of  the  solutions  to  the  scattering  problem  for  the  sphere 
are  in  terms  of  spherical  vector  wavefunctions  which  comprise  a  complete  set  of  solu¬ 
tions  to  the  vector  wave  equation  in  spherical  coordinates.  Details  of  these  solutions 
are  given  in  the  text  by  Stratton  (ref  14).  The  problem  of  scattering  by  a  spherical  shell 
was  first  derived  by  Aden  and  Kerker  in  1951  (ref  23).  These  authors  consider  a  spheri¬ 
cal  scatterer  consisting  of  two  concentric  homogeneous  regions  which  are  in  direct 
contact:  a  spherical  interior  region  covered  with  a  concentric  spherical  shell.  The 
general  solutions  are  given  for  the  case  when  the  interior  region  and  the  concentric 
spherical  shell  have  different  complex  dielectric  constants  and  different  complex 
propagation  constants.  The  study  by  Aden  and  Kerker  was  extended  by  Wait  in  1964  to 
the  case  when  the  scatterer  consists  of  any  number  of  concentric  homogeneous  re¬ 
gions  (ref  24). 

Stimulated  by  interest  in  nuclear  electro-magnetic  pulse  (EMP),  studies  of  the 
shielding  of  transient  fields  by  a  spherical  metalic  shield  were  conducted  in  1965  by 
Harrision  and  Papas  (ref  25).  Using  the  steady  state  solutions  for  scattering  by  an 
incident  planewave  of  frequency  f  (ref  14),  the  authors  calculate  the  field  amplitudes  of 
the  incident  wave,  that  is,  the  shielding  effectiveness,  as  a  function  of  frequency.  This 
result  is  just  the  transfer  function  for  the  field  amplitudes,  either  electric  or  magnetic,  at 
the  center  of  the  sphere  due  to  an  incident  wave  oi  frequency  f.  Using  this  transfer 
function  and  standard  Fourier  transform  techniques,  these  authors  make  numerical 
calculations  of  the  transient  fields  at  the  center  of  the  sphere  due  to  a  Gaussian  EMP 
wave  form  incident  on  the  sphere.  Further  investigations  of  the  steady  state  shielding 
problem  were  made  by  Chui,  Dudley,  and  Bristol  in  1969  (ref  26).  These  authors  use 
the  vector  multipole  expansions  of  the  fields  described  in  Jackson  (ref  27)  instead  of 
Stratton’s  vector  spherical  wavefunctions.  Mathematically,  however,  these  two  formula¬ 
tions  are  equivalent.  The  only  other  difference  in  their  approach  is  that  they  consider  a 
circularly  polarized  incident  wave  rather  than  a  linearly  polarized  wave,  which  adds 
some  symmetry  to  the  problem.  This  work  extends  the  results  of  Harrison  and  Papas 
by  calculating  the  magnitudes  of  the  electric  and  magnetic  fields  in  the  interior  of  the 
spherical  shell  along  the  z-axis,  and  not  just  at  the  center  z=0.  These  calculations  were 
performed  at  six  frequencies,  including  the  first  two  resonant  frequencies.  In  addition, 
they  present  plots  of  the  equi-field  curves  throughout  the  interior  region  for  two  different 
frequencies.  The  results  of  this  work,  and  additional  work  by  Dudley  and  Quintenz  in 
1975  (ref  28)  show  that  the  magnitudes  of  the  fields  away  from  the  center  of  the  sphere 
may  be  four  orders  of  magnitude  greater  than  the  fields  at  the  center. 

The  theory  of  the  shielding  effectiveness  of  a  spherical  conducting  shell  is  devel¬ 
oped  in  this  report  with  reference  to  the  fields  at  the  center  of  the  sphere.  The  approach 
used  here  is  the  same  as  that  of  Harrison  and  Papas  (ref  25).  Numerical  calculations 
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are  performed  for  a  sphere  with  outer  diameter  d=155  mm  and  thickness  b-a=0.1  mm, 
over  a  range  of  frequencies  from  10°  to  1010  Hz.  The  results  show  a  series  of  discrete 
resonances  at  high-frequencies.  Although  these  resonances  are  analogous  to  those  for 
the  cylinder  problem,  the  have  apparently  not  been  thoroughly  investigated  in  previous 
work.  Equations  for  the  resonance  conditions  are  derived.  In  addition,  an  approximate 
formula  for  the  shielding  effectiveness  is  obtained  which  is  accurate  for  frequencies  in 
the  resonance  regime.  To  this  author’s  knowledge  this  result  has  not  been  derived 
previously. 

The  contents  of  this  report  are  organized  as  follows.  The  theory  and  results  for  the 
cylinder  are  presented  in  part  1 ,  and  the  theory  and  results  for  the  sphere  are  presented 
in  part  2.  In  the  first  section,  section  1 ,  the  field  equations  will  be  solved  by  reducing  the 
Maxwell  equations  to  a  scalar  wave  equation,  the  scalar  Helmholtz  equation,  in  cylindri¬ 
cal  coordinates.  In  section  2  and  section  3  the  exact  solutions  for  the  TM  and  TE  fields 
shall  be  presented  in  the  forms  of  series  expansions  of  the  fields  and  the  expansion 
coefficients  determined  by  the  application  of  boundary  conditions  on  the  cylindrical 
surface.  In  sections  4,  5,  and  6  the  field  equations  shall  be  solved  for  the  case  of  an 
obliquely  incident  planewave,  that  is,  a  planewave  incident  at  an  arbitrary  angle  to  the 
z-axis.  Then  in  section  7  the  shielding  effectiveness  shall  be  defined  and  the  results  of 
numerical  calculations  of  the  shielding  effectiveness  versus  frequency  for  the  cylinder 
shall  be  presented  for  both  the  E-polarized  and  H-polarized  incident  waves.  The  nu¬ 
merical  methods  for  evaluation  the  Bessel  functions  are  described  in  Appendix  A.  The 
results  are  compared  to  the  approximate  formula  for  the  shielding  effectiveness  which  is 
derived  in  appendixes  B,  C,  and  D.  An  alternate  definition  of  shielding  effectiveness  is 
defined  in  appendix  E.  In  part  2,  the  Maxwell  equations  are  solved  by  using  the  well 
known  spherical  vector  wavefunctions  which  form  a  complete  set  of  solutions  to  the 
vector  wave  equation,  the  vector  Helmholtz  equation,  in  spherical  coordinates.  In 
section  8  the  scattering  of  a  planewave  by  a  spherical  conducting  shell  is  solved  by 
expanding  the  electromagnetic  fields  in  terms  of  spherical  vector  wavefunctions.  The 
expansion  coefficients  are  obtained  by  matching  the  boundary  conditions  at  the  spheri¬ 
cal  surface.  In  section  9  the  field  solutions  are  used  to  obtain  expressions  for  the 
shielding  effectiveness  of  the  sphere,  and  in  section  10  the  numerical  calculations  of  the 
shielding  effectiveness  are  presented.  The  resonance  conditions  and  the  approximate 
formula  for  the  shielding  effectiveness  are  derived  in  appendix  F.  Finally,  the  conclu¬ 
sions  are  summarized  in  section  1 1 . 

THEORY:  PERPENDICULAR  INCIDENCE 

An  infinitely  long  cylindrical  conducting  shell  of  inner  radius  r=a  and  outer  radius 
r=*b  has  electrical  conductivity  a,  dielectric  constant  e,  and  magnetic  permeability  p.  It 
is  assumed  that  the  conductor  is  a  non-magnetic  material,  that  is,  it  is  a  diamagnetic  or 
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paramagnetic  material,  such  as  copper  or  aluminum,  and  therefore  p=p0.  The  case 

when  the  conductor  is  a  ferromagnetic  material  is  more  complex  due  to  the  frequency 
and  magnetic  field  dependence  of  the  permeability,  and  therefore  shall  not  be  consid¬ 
ered  here.  Consequently,  this  analysis  does  not  rigorously  apply  to  shielding  materials 
such  as  iron  or  steel.  It  shall  be  assumed  in  the  analysis  that  p=p0- 

The  cylinder  is  located  in  empty  space  (vacuum)  which  has  a  dielectric  constant  eo 
and  permeability  po.  A  linearly  polarized  planewave  of  frequency  f  is  incident  on  the 

cylinder  with  its  direction  of  propagation  perpendicular  to  the  cylinder  axis.  Two 
polarization  states  shall  be  analyzed  separately:  one  state  of  linear  polarization  of  the 
incident  wave  with  the  E-vector  parallel  to  the  cylinder  axis  (TM),  and  one  polarization 
state  with  the  H-vector  parallel  to  the  cylinder  axis  (TE).  Any  state  of  elliptical  polariza¬ 
tion  of  the  incident  planewave  may  be  expressed,  in  general,  as  a  sum  of  these  two 
polarization  states,  TM  and  TE,  with  complex  coefficients.  Therefore,  once  solutions  for 
the  TM  and  TE  incident  waves  are  found,  the  general  solution  for  an  arbitrary  state  of 
polarization  of  the  incident  wave  may  be  obtained  by  superposition  of  these  independ¬ 
ent  solutions. 
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Figure  1.  Geometry  for  perpendicular  incidence.  Top  view. 


The  physical  geometry  for  this  problem  is  shown  in  figure  1. 
There  are  three  distinct  regions:  I,  II,  and  III,  for  which  the 
electromagnetic  fields  must  be  determined .  The  boundary 
conditions  on  the  field  vectors  at  the  interfaces  of  the 
different  regions  are  as  follows:  the  tangential  E-fields  and 
tangential  H-fields  must  be  continuous,  and,  the  normal  D-fields 
and  normal  B-fields  must  be  continuous.  The  electromagnetic 
field  equations,  the  Maxwell  equations,  in  MKSA  units  for  a 
general  conducting  medium  are  as  follows: 


VxE 


dB 

3T' 


V*H 


J  + 


dD 

cTE' 


V • D  =  0, 

V«B  =  0, 


where  the  material  equations  are 


D  =  EE, 
B  =  MU, 
J  =  OE. 


Since  the  incident  planewave  has  a  harmonic  time  dependence,  thi 
implies  that  all  the  fields  shall  have  the  same  time  dependence. 


Therefore  it  is  assumed  that  the  field  vectors  E,  H,  D,  and  B 

vary  in  time  as  e_ico^.  Using  this  fact  together  with  the 
constitutive  equations,  the  Maxwell  equations  take  the  form 


VxE  * 

icoy  0H, 

(1) 

VxH  = 

(O-ittE )E, 

(2) 

V*E  * 

0, 

(3) 

V  »H  = 

0. 

(4) 

The  wave  equations  for  the  E  and  H-fields  may  be  derived  from 
the  field  equations.  Taking  the  curl  of  equation  1  and  using  the 
vecto-r  identity 

Vx(VxE)  =  V ( V«E)  -  V *E 

and  then  substituting  equations  2  and  3  into  the  result  yields 
the  vector  wave  equation,  or  Helmholtz  equation. 


where 


( V  a+y a )E  =  0, 
ya  =  icopo  (o-i«E )  . 


(5) 


Similarly,  by  taking  the  curl  of  equation  2,  and  then 
substituting  equations  1  and  4,  the  wave  equation  for  the  H-field 
takes  the  same  form: 


where 


( Va+ya )H  =  0, 
ya  =  icoMo  (o-iwE ) . 


(6) 

(7) 


In  the  case  when  the  conductivity  is  zero,  that  is,  in  regions  I 
and  III,  the  same  wave  equation  holds  with  o= 0,  and  e  =  e0,  which 
implies  y-k.  As  a  mathematical  note  it  should  be  mentioned  that 
the  wave  equations,  equations  5  and  6,  are  automatically 
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satisfied  by  the  E  and  H-fields  if  the  fields  satisfy  Maxwell's 
equations.  That  is.  If  E  and  H  satisfy  Maxwell's  equations  then 
they  satisfy  5  and  6.  However,  if  the  fields  satisfy  5  and  6  this 
does  not  imply  that  they  satisfy  the  Maxwell  equations. 

Therefor©,  once  solutions  of  the  wave  equations  are  obtained  it 
is  necessary  to  verify  that  they  satisfy  the  Maxwell  equations. 

In  the  first  case  to  be  considered  the  incident  E-field  is 
polarized  in  the  z-direction: 

E*nc  -  E„eikx. 

Symmetry  about  the  z-axis  implies  that  the  total  E-field  in  each 
of  the  three  regions  is  also  polarized  in  the  z-direction.  Thus 
the  only  nonzero  component  of  the  E-field  is  Ez  in  the  TM  case. 

By  using  the  first  Maxwell  equation,  equation  1,  and  the 
expression  for  the  curl  operator  in  cylindrical  coordinates,  the 
total  H-field  has  two  nonzero  components  given  by 

1  1  dEz 

Hr  ISpT  FoRT'  <8> 


Thus,  when  the  electric  field  component  Ez  is  known  the  magnetic 

fields  may  be  obtained  by  differentiation.  The  electric  field  Ez 

will  be  obtained  in  each  of  the  three  regions  I,  II,  and  III  by 
solving  the  respective  wave  equations.  The  magnetic  fields  will 
be  calculated  from  Ez  by  using  equations  8  and. 9.  By  using  this 

approach  Maxwell’s  equations  are  automatically  satisfied. 
Similarly,  in  the  TE  case  the  only  nonzero  component  of  the  H- 
field  is  Hz .  Using  the  second  Maxwell  equation,  equation  2,  the 

components  of  the  E-field  in  this  case  are  given  by 

1  1  dHz 

Er  "  r  cft“'  (10) 

-1  dHz 

E<t>  *  37"'  (11) 
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The  material  parameters  in  region  I  are  y0»  and  o=0.  The 

wave  equation  for  E  (TM  case)  or  H  (TE  case)  is  given  by 

z  z 


( V*+k*  )\J>  «  0, 


k* 


. 

c* ' 


where  \)>  is  equal  to  E  or  H  .  Using  the  expression  for  the 

Z  Z  ^ 

Laplacian  in  cylindrical  coordinates  this  becomes 


1  d 
r  Hr 


+  r*H^  + 


k*\|> 


-  0  . 


(12) 


To  solve  this  equation  use  separation  of  variables.  Let 
\|>=R(r)Q(<t>) .  Then  the  equation  for  the  angle  variable  becomes 


aQ 

<T<F 


+ 


n*Q 


0, 


(13) 


which  has  solution 


Q(<t>) 


«±ln*. 


The  parameter  na  is  the  separation  constant.  To  ensure  that  the 
solution  is  single  valued  n  must  be  an  integer:  n=  0,±1,±2,... 
The  radial  equation  takes  the  form 


d*R  ,  1  dR  , 

Hr7  r  3r  + 


(k2+  pi) R 


0. 


(14) 


If  the  substitution  x=kr  is  made,  then  this  equation  takes  the 
standard  form  of  Bessel's  differential  equation: 


daR 
dx 2 


+ 


dR 

3x 


+  (1  - 


n* 

x1 


•)R 


o. 


(15) 


The  linearly  independent  solutions  are  the  Bessel  functions  of 
the  first  and  second  kinds 


R(x)  =  Jn(x),  Yn(x). 
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Another  set  of  linearly  independent  solutions  are  the  Hankel 
functions 

HA(X)  "  Jn(x)  +  iYn(x)' 

HA(X)  "  Jn(x)  "  iYn(x)* 

The  total  field  in  region  I  may  be  written  as  a  superposition  of 
the  incident  and  scattered  fields  in  the  form 


*  «ikX  +  l  in  V^tkr) 


e 


in<t> 


(16) 


n=— 


The  Hankel  function  H‘(kr)  is  chosen  since  it  represents  an 
outgoing  wave  for  large  r.  The  incident  field  may  be  expanded  in 
a  Fourier  series  as  follows 


e 


ikx  „ikrcos(<J>) 


£  in  Jn(kr)  ein4> 


n»- « 


(17) 


Therefore,  using  this  expansion  the  total  field  may  be  written 


y  =  l  in[Jn(kr)  +  anH^(kr )  ]ein<^,  r>a  . 


(18) 


n=-« 


Next,  solutions  shall  be  found  for  the  wave  equation  in  region 

II.  The  material  parameters  in  region  II  are  y0,  e,  and  o. 

The  wave  equation  for  E  in  the  TM  case,  or  H_  in  the  TE  case,  is 

z  z 


where,  from  7 


(19) 
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As  in  the  previous  derivation,  separation  of  variables  leads  to 
the  solution 


♦  -  l  in  [bnJn(yr)  +  cnHA(yr)^in<l>  '  a<r<b  •  (20) 

n» 


The  combination  of  functions  Jn(yr)  and  H^(yr)  are  chosen  so  the 

solution  in  this  region  has  a  similar  form  to  the  solution  in 
region  I.  The  arguments  of  these  functions  are  now  complex, 
however,  since  y  is  complex.  For  a  good  conductor,  and  for 
frequencies  such  that  {o/we0)>>l,  approximate  relations  for  y 
are,  from  19, 


§;  (i 


(OE  o 


and 


‘60 


In  region  III  the  material  parameters  are  vi0,  e0,  and  o=0. 
Thus  the  wave  equation  has  the  same  form  as  in  region  I: 

( V*+k2 ) \p  -  0,  k2  -  |I. 


As  before,  separation  of  variables  yields  the  solution 


¥  -  5  in  dnJn(kr)eln<t>'  r<a.  (22) 


Only  the  Bessel  function  Jn(kr)  is  used  since  the  other 

solutions,  which  contain  Neumann  functions,  approach  infinity  as 
r  approaches  zero.  These  other  solutions  are  physically 
unacceptable. 

Due  to  the  fact  that  only  the  Hankel  function  H^(x)  shall  be 

used,  the  superscript  (l)  shall  be  omitted  hereafter.  Therefore, 
whenever  the  function  H  (x)  appears  it  should  be  understood  that 

it  stands  for  H*(x).  The  reason  for  this  is  to  avoid  confusion 
n 

with  the  prime  symbol  ( ' )  which  is  used  for  differentiation.  This 
convention  shall  be  adhered  to  throughout  this  thesis. 
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2.  TM  INCIDENT  WAVE 

The  full  solutions  for  the  fields  in  the  TM  case  are  obtained 
by  using  equations  20,  21,  and  22  which  give  solutions  for  E  and 

Z 

then  using  equations  8  and  9  to  determine  Hr  and  H^.  The  full 
electromagnetic  field  solutions  are,  in  region  I, 


E2  -  E0  l  in[Jn(kr)  +  anHn(kr)]ein<t>, 


r>b; 


n«— 


E0  1 


H_  *  -T^  ±  l  nin+1[Jn(kr)  +  anHn(kr)]eintp,  r>b; 

n=-« 


,  in<t» 


r  icojj  0 


H<t>  *  SS?  S  intJA<kr>  +  anHA<kr>3e 

n=-« 

region  II: 


in<|> 


r>b; 


Ez  -  E0  l  in[bnJn(yr)  +  cnHn(yr >  ' 

n=-<*> 


a<r<b; 


H. 


r  -  lifr;  i  J  nin+1[bnJn(Yr)  +  =„Hn(yr) ]eln®.  a<r<b; 


,n+i 


in<f> 


n=-<*> 


in[VA^r>  +  cnHA(Yf)]eln,,>'  a<r<b-' 


n=-» 


region  III: 


Eo  i  in  dnJn(kr)ein<t>' 
n=-«» 


r<a; 


Hr  T§t7  ?  S  ni"+1  anJn<kr>' 
n= 


,in<t> 


r<a; 


H 


-kE, 


l  in  <VA{kr)ein<t>' 


4>  icon , 


r<a , 


n=-«° 


(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 
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The  coefficients  an,  bn,  cfl,  and  dn  are  determined  by  matching 
boundary  conditions  at  the  interfaces.  Continuity  of  the 


tangential  E-field  at  r-b  implies 

Jn(kb)  +  anHn(kb)  "  bnJn(yb)  +  cnHn(yb)*  (32) 

Continuity  of  the  tangential  E-field  at  r-a  implies 

br,J«(Ya)  +  CnH«(Ya)  “  dr,Jr,<ka)-  (33) 

n  n  '  n  n  '  n  n  ' 

Continuity  of  the  tangential  H- field  at  r«b  implies 

k[j^(kb)  +  anH^(kb)]  =  Y[bnJ^(yb)  +  cnH^(yb)].  (34) 

Continuity  of  the  tangential  H-field  at  r*a  implies 

y[bnJn(ya)  +  cnHA(ya)]  "  kdnJn(ka)*  (35) 


These  are  all  the  required  boundary  conditions.  Since  y=y0,  the 
boundary  condition  on  the  normal  component  of  B  is  redundant. 
These  equations;  32,  33,  34,  and  35;  comprise  a  nonhomogeneous 
system  of  four  equations  in  the  four  unknowns  an,  bn,  cfl,  and  dfl. 

The  solutions  of  this  system  are  obtained  directly  with  the  use 
of  Kramer's  rule.  This  results  are  given  by  the  relations 


(36) 


where  the  and  V  are  determinants  given  as  follows: 
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Hn(kb)  -Jn(Yb)  -Hn(Yb)  -Jn(kb> 

0  Jn(ya)  Hn(Ya)  0 

U4  -  ,  (40) 

kH^(kb)  -yJ^(yb)  -yH^(yb)  -kJ^(kb) 

0  YJn^ya)  YH^(ya)  0 


and 


Hn(kb)  -Jn(Yb)  -Hn(yb)  0 

0  Jn(ya)  Hn<Ya)  -Jn(ka) 

v  -  .  (41) 

kH^(kb)  -yJ^(yb)  -yH^yb)  0 

0  YJ^ya)  yH^(ya)  -kJ^(ka) 


This  completes  the  solutions  for  the  fields  in' the  TM  case.  The 
expansion  coefficients  may  be  determined  explicitly  by  evaluating 
the  above  determinants. 

Although  the  series  solutions  of  the  fields  are  rather 
cumbersome  since  the  expansion  coefficients  are  given  by 
quotients  of  determinants  of  Bessel  and  Hankel  functions,  they  do 
constitute  exact  solutions,  and  they  are  sufficient  for  the 
purposes  of  this  study.  It  is  very  desireable,  and  it  seems  very 
reasonable  to  expect  that  the  field  solutions  may  be  expressible 
in  a  closed  form  in  terms  of  an  integral.  This  author  is  not 
aware  if  such  a  closed  form  integral  solution  exists  or  has  been 
deri/ed  previously.  This  is  a  point  of  departure  for  further 
investigations . 
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3.  TE  INCIDENT  WAVE 

The  full  solutions  for  the  fields  in  the  TE  case  are  obtained 
by  using  equations  20,  21,  and  22  which  give  solutions  for  H  and 

Z 

then  using  equations  10  and  11  to  determine  Er  and  E^.  The  full 
electromagnetic  field  solutions  are,  in  region  I: 


Hz  -  H0  J  in[Jn(kr)  +  AnHn(kr)]ein4>, 


r>b; 


r  -  ISeT  F  s  nln,'1[Jn(kr)  +  AnHn(kr)  ]eln<p,  r>b; 


n= 


4>  -  f£tr  intJA<kr)  +  ^(krlje1"®.  r>b; 


n=-o» 


region  II: 


oo 

Hz  -  H„  J  in[BnJn(yr)  +  CnHn(yr)]ein't',  a<r<b;  (44) 


Er  '  X  nin+1[BnJn(yr)  +  CnHn(yr)  ]eln,,,,  a<r<b;  (45) 

n=-«*> 


oo 

E4,  -  irtfe  E  ‘n[v;(yr)  +  ^(YDje^,  a < r < b ;  (46) 


n*-« 


region  III: 


Hz  -  H 0  l  in  DnJn(kr)ein<p, 
n--« 


r<a; 


-Ho  1 


'r  i«E, 


i  l  ni  DnJn(kr)«"*v, 


r<a; 


'<t>  “  iwe, 


l  in  DnJ^(kr)ein<t>, 


r<a . 


n= -oo 
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As  in  the  TM  case,  coefficients  An,  Bn,  Cn,  and  Dn  are 

determined  by  matching  boundary  conditions  at  the  interfaces. 
Continuity  of  the  tangential  H-field  at  r-b  implies 


Jn(kb)  +  AnHn(kb)  *  BnJn(yb)  +  cnHn(yb)*  (50) 


Continuity  of  the  tangential  H-field  at  r=a  implies 


B  J { y a )  +  CH  (ya)  -  D  J  (ka). 
n  n  '  nnw  nn 


Continuity  of  the  tangential  E-field  at  r=b  implies 


TSz:  tJA(kb)  +  AnHA<kb>]  *  tBnJA(yb)  +  CnHA<yb>] 


Continuity  of  the  tangential  E-field  at  r-a  implies 


O—  1  WE 


[VA‘V«>  +  CnHA<^>]  -  T§F7  DnJA(ka) 


These  are  all  the  required  boundary  conditions.  These  equations; 
50,  51,  52,  and  53;  comprise  a  nonhomogeneous  system  of  four 
equations  in  the  four  unknowns  An,  Bn,  Cn,  and  Dn<  As  before, 

the  solutions  are  obtained  with  the  use  of  Kramer's  rule.  The 
results  are  as  follows: 


where  the  X.^  and  Y  are  determinants  given  as  follows: 
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0 


X, 


X, 


X, 


-Jn(kb) 


ISF7JA<kb>  cdhJA<Yb>  5rfeeHA(Yb)  o 


-Jn<yb) 


Jn(Y») 


-Hn(yb) 


Hn(ya) 


-Jn(ka) 


^fejA<va>  5^eha<y») 


(55) 


Hn(kb) 


-Jn(kb) 


-k 


-Hn(yb) 


Hn(ya) 


-Jn(ka) 


I^7HA<kb)  TSi7JA(kb)  o^feHA<vb>  0 

I^7JA(ka) 


o-iwe  n 


(56) 


Hn(kb) 


-Jn(Yb) 


Jn<ya) 


-Jn(kb) 


ISF7HA>kb'  dfeJA<Yb)  T^F7JA(l'b)  0 

X - 7  'Ival  0  ,  ~"J  '  (ka) 


-Jn(ka) 


o-iweJn^a^ 


iwe  -  n 


(57) 
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This  completes  the  solutions  for  the  fields  in  the  TE  case.  The 
expansion  coefficients  may  be  determined  explicitly  by  evaluating 
the  above  determinants. 

This  concludes  the  presentation  of  the  electromagnetic  field 
solutions  for  the  scattering  of  a  normally  incident  planewave 
from  a  circular  cylindrical  conducting  shell  of  inner  radius  r=a 
and  outer  radius  r*b.  The  expansion  coefficients  of  the  fields  in 
the  region  r<a  shall  be  evaluated  in  the  Appendixes  for  the 
purposes  of  calculating  the  shielding  effectiveness.  In  the  next 
sections,  sections  4,  5,  and  6  the  solutions  for  scattering  by  an 
obliquely  inciwjnt  planewave  shall  be  derived. 
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4.  OBLIQUE  INCIDENCE 


In  general.  Maxwell's  equations  may  be  solved  by  the  use  of 
two  auxilary  vectors  n,  and  IT,,  called  Hertz  vectors.  The  time 
harmonic  Maxwell  equations,  equations  1,  2,  3,  and  4,  are 
automatically  satisfied  if 


(60) 

(61) 

and  if  the  Hertz  vectors  satisfy  the  vector  wave  equations 


E  =  VxVxrij  +  iwy  VxTIj 
H  ■  (o-iwEjVxlTj  +  VxVxir,, 


vxvxn,  -  v(v«n,)  -  y*ni  -  o, 
VxVxna  -  v(v*na)  -  y*na  •  o, 

where,  as  before, 

ya  *  iwp0 (o-iwe ) . 


(62) 

(63) 


(64) 


This  may  be  verified  by  inserting  equations  60  and  61  into 
equations  1  and  2,  and  then  using  equations  62  and  63.  In  the 
case  of  cylindrical  symmetry  the  only  nonzero  components  of  the 
Hertz  vectors  are  the  z-components.  If  n,*¥,e  and  n*=0,  where 

is  a  scalar  function,  then  the  cooresponding  solutions  of 
Maxwell's  equations  are  given  by 


E  =  VxVxn,  =  ih  er  +  ih  i  e0  +  (ya-  h*)T,  ez,  (65) 

H  -  (o-i(OE)Vxn,  *  (o-icoE)  (■p  er  -  e^),  (66) 


<  Jl2 

where  it  is  assumed  that  -  \|>,(r,<J>)e  .  A  z-dependence  of  the 

jkhz 

form  e  is  chosen  to  match  the  form  of  the  incident  wave.  The 
vectors  er,  e^,  and  e2  are  the  usual  unit  vectors  in  the 
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cylindrical  coordinate  system.  An  independent  set  of  solutions  of 

Maxwell's  equations  is  obtained  if  n,=0  and  n2=¥2e  which  yields 

£* 

E  =  icoy  Vxn,  =  iwy  (i  er  -  ^  e^),  (67) 

H  -  VxVxn,  =  ih  §r  er  +  ih  r  Eb  e4>  +  (y*"  h’>*2  ez'  (68) 


4 

where  ¥2  =  \|>2  (r,4>)e  .  Note  that  the  solutions  based  on  the 

wavefunction  S',  have  a  nonzero  E  component,  and  the  solutions 

Ct 

based  on  the  wavefunction  ¥2  have  a  nonzero  H2  component. 

Therefore,  equations  65  and  66  represent  TM  fields  and  equations 
67  and  68  represent  TE  fields. 

For  the  case  of  perpendicular  incidence  a  TM  incident  wave 
gives  rise  to  a  TM  scattered  wave  and  a  TE  incident  wave  gives 
rise  to  a  TE  scattered  wave.  This  is  not  true  in  the  case  of 
oblique  incidence.  For  oblique  incidence  a  pure  TM  incident  wave 
gives  rise  to  a  combination  of  TM  and  TE  scattered  vaves.  And, 
similarly,  a  pure  TE  incident  wave  gives  rise  to  a  combination  of 
TE  and  TM  scattered  waves.  Therefore  both  TM  and  TE  components  of 
the  scattered  fields  are  required  to  express  the  total  fields  in 
each  of  the  regions  I,  II,  and  III.  This  implies  that  both 
wavefunctions  ¥,  and  must  be  used  to  give  independent  TM  and 
TE  solutions  for  the  scattered  fields,  respectively. 

In  region  I,  r>b,  the  material  parameters  are  z0,  y0,  and  o=0. 
The  scalar  Helmholtz  equation  takes  the  form 


(V2  +  k2  )¥  =0,  k2  =  (69) 


The  solutions  for  the  wavefunctions  are 


(TM) 

-  |f  eihz  J  in  anHn(E,r)ein<t>, 

-oo 

(70) 

(TE) 

'Hz  =  eihz  £  in  bnHn(Jlr)ein<t), 

—  oo 

(71) 
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2  2  2 

where  8,  «k  -  h  .  The  Hankel  function  is  chosen  so  that  the 
scattered  field  reduces  to  an  outgoing  wave  as  r  approaches 
infinity.  The  constants  E0/kH  and  H„/ki  are  included  so  that  the 
expansion  coefficients  are  dimensionless.  For  convenience,  the 
field  strengths  E0  and  H„  will  be  chosen  to  be  equal  to  those  of 
the  incident  planewave.  In  region  II,  a<r<b,  the  region  within 
the  conductor,  the  material  parameters  are  e,  y0,  and  o.  The 
Helmholtz  equation  takes  the  form 


( V*+y * )T  -  0,  y*  -  iwyc(o-i«e). 


The  solutions  for  the  wavefunctions  in  region  II  are 


(TM) 

*i  -  h  «lhZ  5  in[cnJn<*r>  +  dnHn(xr)]ein<t>' 

(72) 

(TE) 

*.  -  jjf  eihZ  5  lntVr.<V'r|  +  9nHn<xr)]el  t 

^  —CO 

(73) 

2  2  2 

where  x  =y  -  h  .  Both  Bessel  and  Hankel  functions  are  required  in 
the  expansions  72  and  73  in  order  to  match  the  interior  and 
exterior  fields.  In  region  III,  r<a,  the  material  parameters  are 
e„  and  p0.  Therefore,  the  wavefunctions  in  the  interior  region, 
region  III,  take  the  forms 


(TM)  *.  «  §f  eihz  2  in  PnJn(fcr)ein<t>‘,  (74) 

(TE)  -  |f  eihz  2  in  qnJn(fcr)ein<l>.  (75) 

—  CD 


Only  the  Bessel  function  is  required  since  it  is  finite  at  r=0. 
Using  these  results  for  the  wavefunctions  the  scattered  fields  in 
a  given  region  may  be  obtained  by  using  equations  65,  66,  67,  and 
68.  For  example,  in  the  exterior  region,  region  I,  where  r>b,  the 
TM  scattered  fields  art  obtained  by  inserting  the  wavefunction  70 
into  the  set  of  equations  65  and  66,  and  the  TE  scattered  fields 
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are  obtained  by  inserting  71  into  equations  67  and  68.  The 
resulting  TM  scattered  fields  are,  for  r>b. 


And, 


E™  =  iE0eihz  l  in  £  anH^(8-r)ein<t>, 

—co 

„TM  ihz  n  ,  n  h  , ,  _  .  =  „  ,  „in<t> 

E<j>  =  iE0e  l  i  ^(in)  anHn(fcr)e 

— 09 

E™  «  E0eihz  i  in  I  anHn(8,r)ein4>' 

-09 

H™  =  -iH0eihz  5  in  anHn(ar)ein<t>' 

—  oo 

H™  =  iH0eihz  l  in  anH^(4r)ein<,>. 

—  oo 

the  resulting  TE  scattered  fields  in  the  region  r>b  are 
Ef  =  iE0eihz  l  in  ||  bnHn(J'r)ein<,>' 

—  oo 

E™  .  -iE„elhz  i  i"  bnH’(lr)ein*, 

—  OO 

h£E  =  iH0eihz  l  in  £  bnH^(ar)ein<t>, 

— oo 

hJE  -  iH0eihz  l  in  bnHn(ar)ein<t>, 

— oo 

UTE  u  ihz  n  .n  %  K  u  .  n  .  in<t> 

Hz  =  H0e  2j  1  £  bnHn(J,r)e 

—  oo 
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The  scattered  fields  in  the  other  regions  are  obtained  similarly. 
In  general,  the  scattered  fields  are  of  the  form 


E 


s 


,TM 


+  E 


,TE 


Having  determined  explicit  expressions  for  the  scattered 
fields  the  next  step  is  to  obtain  expressions  for  the  incident 
fields.  The  incident  planewave  propagates  at  an  angle  0  to  the  x- 
axis  as  shown  in  Figure  2.  The  k-vector,  or  wave  normal,  lies  in 
the  x— 2  plane.  Two  distinct  states  of  linear  polarization  of  the 
incident  wave  will  be  treated  seperately.  The  case  in  which  the 
E-vector  of  the  incident  wave  lies  in  the  x-z  plane  shall  be 
called  the  TM  incident  wave.  And  the  case  in  which  the  H-vector 
of  the  incident  wave  lies  in  the  x-z  plane  shall  be  called  the  TE 
incident  wave.  In  rectangular  coordinates  the  electromagnetic 
fields  of  the  TM  incident  wave  are  given  by 


(TM) 

Einc 

(TM) 

Hinc  ! 

And  the 

fields  of 

(TE) 

Einc 

(TE) 

Hinc 

where  H0 

e  o / y  o  E 

_  .  ^  ^  ~  _  .  ik(xcos0  +  zsin0) 

*  Eo(-sm0  e  +  cos0  e  )  e  '  , 

X  z 


-H„e  e 


ik(xcos6  +  zsin0) 


E0e..  e 


ik(xcose  +  zsin0) 


inc  =  H0 ( -sine  ex  +  cos©  e2)  e 


ik(xcos0  +  zsin©) 


wave  in  cylindrical  coordinates  may  be  determined  by  using  the 
transformation  of  basis  vectors: 


e  =  cos<t>  e  +  sin0  e 
r  x  y 


e^  -sin0  ex+  cos(|)  e^. 
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Figure  2.  Geometry  for  oblique  incidence.  Side  view. 


The  cylindrical  components  of  the  TM  incident  wave  are 


(TM) 


E.(-sinecos<Melhzeikxcose, 


_  . „ihz  ikxcosG 

E0  (sin6sin<t>)e  e 


„  / \  _ihz_i3cxcos0 
E0(cos8)e  e 


Hr  = 


„  ,  „ihz  ikxcosG 

H0(-sino)e  e 


H»  ■ 


u  /  . „ihz  ikxcosG 

H0  ( -cos4> )  e  e 


where  h-ksinG.  Similarly,  the  cylindrical  components  of  the  TE 
incident  wave  are 


(TE) 


it  o  \  — ihz  ikxcosG 

Er  ■  E0(sino)e  e 


_  _  . _ „ . ,  ihz  ikxcosG 

E^  =  E0  ( cos<t> )  e  e 
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Hr  -  H0(-sinecos4»eihzeikxcose( 


-  H0  (sin9sin<t>)eihzeikxcos0, 


H2  «  H0 (cose )eihzeikxcos0 


The  incident  fields  may  be  expanded  in  a  series  of  Bessel 
functions  by  using  the  Fourier  series  representation 


,ikrcos<t>cos0 


l  in  Jn(fcr)ein4), 


where  8,=kcos0 .  The  following  identities  are  also  needed 

,  .ikrcos<J>cos0  1  d  /„ikrcos<t>cos0  i 

sin<t>  e  -  -n?  3*  >' 

cos<t>  ei^rcos4)cose  =  (.i)l  d_  (eikrcos<i>cos0  j  # 

where  fi,=kcos0.  Using  these  identities  the  required  expansions  of 
the  TM  incident  wave  are  as  follows 


(TM) 

Er  " 

E„elhz 

E  in  (i>£  J^(&r)ein4>, 

— *o 

(76) 

(TM) 

E*  - 

E0eihz 

S  in  Jn(^r)ein<t>, 

— oo 

(77) 

(TM) 

Ez  " 

E.elhz 

5  inf  Jn(&r)ein<t>, 

(78) 
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(TM) 

Hr  * 

H0eihz  £  in 

-  oo 

(79) 

(TM) 

H4»  " 

H„elhz  jj  in 

(i)  Jn(Hr)ein<t>. 

(80) 

Similarly, 

the  required  expansions  of  the  TE  incident  wave 

are 

(TE) 

Er  " 

-E0eihz  £  in 

— oo 

Hr «V*r)ein<l>, 

(81) 

(TE) 

E*  " 

-E0eihz  l  in 

—  oo 

(i)  Jn(ar)ein<>, 

(82) 

(TE) 

Hr  " 

H0eihz  l  in 

— oo 

(i)£  J^(ar)ein<t>, 

(83) 

(TE) 

H*  " 

H0eihz  l  in 

— oo 

(i)icnr(in)  Jn(ftr>ein4>' 

(84) 

(TE) 

Hz  “ 

t t  - ihz  n  j  n 
H0e  i  i 

— oo 

S  Jn(ar>ein4>* 

(85) 

Having  solved  Maxwell's  equations  for  the  scattered  fields 

and 

then  obtained  series  expansions  of  the  incident  planewave,  the 
total  fields  may  be  written  down  at  once.  The  total 
electromagnetic  fields  in  the  region  exterior  to  the  cylinder, 
for  r>b,  are  given  by  the  sum  of  the  incident  and  scattered 
fields.  In  the  interior  region,  for  r<a,  and  inside  the 
conducting  layer,  a<r<b,  the  total  fields  are  given  by  the 
expansions  for  the  scattered  fields  alone.  The  expansion 
coefficients  may  then  be  determined  by  imposing  boundary 
conditions  on  the  fields  at  the  surfaces  of  the  conducting 
cylinder.  This  procedure  shall  be  carried  out  for  both  the  TM  and 
TE  incident  waves  in  the  following  sections,  sections  5  and  6. 
Although  the  full  solutions  shall  be  presented,  the  shielding 
effectiveness  shall  not  be  cax-ulateC  for  the  case  of  oblique 
incidence  due  to  the  complexity  of  the  solutions. 
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5.  OBLIQUE  TM  INCIDENT  WAVE 


In  this  section  the  complete  electromagnetic  field  solutions 
shall  be  obtained  for  the  scattering  of  an  obliquely  incident  TM 
planewave  by  a  circular  cylindrical  conducting  shield.  These 
solutions  differ  from  the  solutions  in  the  case  of  a 
perpendicularly  incident  TM  wave  in  that  the  scattered  fields 
contain  both  TM  and  TE  components.  Therefore,  the  solutions  for 
the  case  of  an  obliquely  incident  planewave  require  eight 
independent  sets  of  expansion  coefficients  instead  of  four.  This 
is  a  consequence  of  the  breaking  of  the  perfect  two-dimensional 
symmetry  in  the  previous  problem.  For  oblique  incidence,  it  is 
interesting  to  note  that  if  the  cylinder  were  a  perfect  conductor 
then  no  mixing  of  the  TM  and  TE  scattered  fields  would  occur  and 
a  TM  incident  wave  would  give  rise  to  a  TM  scattered  wave  only, 
as  in  the  case  of  perpendicular  incidence.  This  implies  that  for 
a  good  conductor  the  secondary  scattered  modes,  TE  in  this  case, 
are  of  low  intensity  relative  to  the  dominant  TM  mode. 

Proceeding  now  with  the  solution.  In  the  exterior  region, 
region  I,  where  r>b,  the  total  fields  are  given  by 


H 


<t> 


H. 


iE.eihz  l  in(£  J^(Hr)  +  £  anH'(lr)  +  ^  bnHn(8,r)  )ein4>,  (86) 

iE0eihz  l  in(^jn(nr)  +  anHn(&r)  -  b^Ur)  )ein<t),  (87) 

—  oo 

E0eihz  l  in(|  Jn(*r)  +  |  anHn(fi,r))ein<t),  .  (88) 

—  oo 

-iH0eihz  £  in(i£jn(s,r)  +  ££anHn(ar)  -  £  bnH^( ar)  )ein4>,  (89) 

— oo 

iH0eihz  l  in( J^( fcr )  +  anH^(&r)  +  R|?( in)bnHn( )  )ein(t> ,  (90) 

—  oo 

H0eihz  2  bnHn(8,rUein<t>-  (91) 
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The  total  fields  inside  the  conductor,  that  is,  in  region  II, 


where  a<r<b,  are 

Er  -  iE6eihz  l  in  {  |  [cnJ^(xr)  +  dnH^(Xr)]  + 

— m 

*  xc  [fnJn(xr)  +  9nHn(’‘r)])eln't’  •  <92> 

E*  ■  iE„elhz  l  in  (]^<in)[cnJn(xr)  +  dnHn(xr)]  + 

-  [fnJn(xr)  +  gnHA(xr)^ein<l>  '  (93) 

E2  =  E0eihz  l  in  £  [cnJn(xr>  +  dnHn(xr)]ein<l>  ,  (94) 

—  GO 

Hr  -  iH0eihz  j  in  [cnJn(xr)  +  dnHn(xr)]  + 

+  £  tfnJA(xr)  +  gnHA(xr)^ein<t>  '  (95) 

H<t>  ■  -iH»eihZ  j  in  UtsI!)  tVA^XD  -  WXD]  ♦ 

-  ^F(in)[fnJn(xr)  +  9nHn(xr)]]ein\  (96) 

Hz  "  H0eihz  l  in  £  [fnJn(xr)  +  9nHn(xr)]ein<t>.  (97) 

—do 


Note  that  the  term  (o-iwe)/(iooE0 )  i*-  equations  95  and  96  is  equal 
to  -y*/k* . 
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The  total  fields  on  the  interior  of  the  shield,  in  region  III, 
where  r<a,  are  given  by 


Er  -  iE,eihz  l  in  (£  PnJ^(Hr)  +  qnJn(Hr))ein*,  (98) 

E  -  iE.eihz  l  in  (^r(in)  PnJn(«-r)  -  qnJ^( Hr )  )ein<t),  (99) 

Ez  -  E0eihz  l  in  £  VnJnll‘r)ein*'  (100) 

— m 

Hr  -  -iH0eihzJ  in  PnJn(S-r)  -  £  qnJ^(Hr))ein<,),  (101) 

%  -  iH0elhz  l  in  (pnJ^(Hr)  +  ^(in)  qnJn( Hr )  )ein4>,  (102) 

Hz  -  H0eihz  l  in  £  qnJn( Hr )eln<^.  (103) 

—  CO 


To  determine  the  expansion  coefficients  the  boundary  conditions 
must  be  applied  at  the  surfaces  r-a  and  r«b.  The  boundary 
conditions  require  that  the  tangential  components  of  the  E-field 
and  the  H-field  be  continuous  on  the  boundary  surfaces.  The 
condition  that  E^  is  continuous  at  r-b  implies 


,^(in)JnUb)  *  ^UniaAl.b)  -  bnH’(lb)  - 

-  in) icnJn(*b)  +  dnHn(xb)J  "  tfnJ’(xb)  +  gnHA(xb)^‘  (104) 


30 


The  condition  that  E  is  continuous  at  r«b  implies 

Z 


|  Vlb>  +  E  «nHn<“»  ■  £  tcnJn‘*b)  +  dnV*b>l-  <105> 

The  requirement  that  is  continous  at  r«b  implies 
J^(lb)  +  anH^(lb)  +  HE(in)bnHn(tb)  - 

£[cnJA<*b>  +  dnHA'Xb)]  ♦  ^(inllVnCXb)  *  3„Hn(xb)].  (106) 


The  requirement  that  Hz  is  continous  at  r*=b  implies 


E  bnHn“b>  -  £  tVn<*b>  +  9nV*b>l 


(107) 


The  requirement  that  is  continous  at  r«a  implies 

E£s(in>tcnV*a>  +  dnV*a>]  “  [fnJA<*a>  +  SnHA<X»>] 


EiS<in>  PnJn(lla)  -  VA(ia> 


(108) 


The  requirement  that  Ez  is  continous  at  r-a  implies 


i  tcnJn<*a>  +  dA'»a|)  ’  E  PnJn * 8,9 ^ 


(109) 


The  requirement  that  HA  is  continuous  at  r=a  implies 
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(110) 


£tcnJA<xa)  +  <SnH^(Xa)]  +  K75a«)[Vn(Xa)  +  %Hn<*a>]  * 

■  PnJi<la>  +  EI5<in> 

And,  the  requirement  that  H  is  continuous  at  r-a  implies 

J  [fnJn(xa)  +  9nHn(xa)]  -  £  Vn|la|'  <1U> 


These  are  all  of  the  required  boundary  conditions.  The  full  set 
of  boundary  conditions  constitutes  a  nonhomogeneous  set  of  eight 
equations  in  the  eioht  unknowns  an,  bn,  cn,  dR,  f  ,  g  ,  p  ,  and 

qfi.  This  system  may  be  solved  by  using  the  standard  technique. 

The  explicit  solutions  for  the  coefficients  is  straightforward 
and  shall  net  be  given  here.  According  to  Kramer's  rule,  the 
solutions  are  given  by  relations  of  the  form 


(112) 


where  the  determinants  and  V  are  determined  by  inspection  of 
the  eight  equations  104,  105,  106,  107,  108,  lC)9,  110,  and  111. 

This  completes  the  solutions  for  the  electromagnetic  fields 
for  the  scattering  of  an  obliquely  incident  TM  planewave.  Due  to 
the  fact  that  the  expansion  coefficients  are  given  by  eight-by¬ 
eight  determinants  the  shielding  effectiveness  for  the  case  of 
oblique  incidence  shall  not  be  calculated. 
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6.  OBLIQUE  TE  INCIDENT  WAVE 


In  this  section  the  complete  electromagnetic  field  solutions 
shall  be  obtained  for  the  scattering  of  an  obliquely  incident  TE 
planewave  by  a  circular  cylindrical  conducting  shield.  As  was 
noted  in  the  last  section  these  solutions  differ  from  the 
solutions  in  the  case  of  a  perpendicularly  incident  TM  wave  in 
that  the  scattered  fields  contain  both  TM  and  TE  components. 
However,  if  the  cylinder  were  a  perfect  conductor  then  no  mixing 
of  the  TM  and  TE  scattered  fields  would  occur  and  a  TE  incident 
wave  would  give  rise  to  a  TE  scattered  wave  only,  as  in  the  case 
of  perpendicular  incidence.  This  implies  that  for  a  good 
conductor  the  secondary  scattered  modes,  TM  in  this  case,  are  of 
low  intensity  relative  to  the  dominant  TE  mode. 

The  expressions  for  the  scattered  fields  in  each  of  the  three 
regions  are  the  same  as  in  section  5  with  the  expansion 
coefficients  replaced  with  capital  letters:  A,  B,  C,  D,  F,  G,  P, 
and  Q.  Only  the  incident  fields  are  different.  The  development  of 
the  solutions  is  identical  to  that  in  section  5.  Proceeding  now 
with  the  solution,  in  the  exterior  region,  region  I,  where  r>b, 
the  total  fields  are  given  by 


Er  -  iE0eihz  l  in(f£  Var>  +  £  AnHA(2,r)  +  t§  bnHn(5,r)  )eluv,  (113) 


E^  =  iE0eihz  l  in(-Jn(S,r)  +  AnHn(&r)  -  Hr )  )einq\  (114) 


E0e  l  i  £  AnHn(S,r)e 


(115) 


Hr  -  -iH0eihz  l  in(-£JA<!'r>  +  I7AnHn(fi,r)  "  K  BnHA( Hr )  )ein4>,  (116) 


4>  “  iH‘eihZ  l  in(£l?Jn<*r>  +  AnHA<*r)  +  KTF  BnHn(l!,r )  (117) 


H0eihz  l  in(jjf  Jn  ( 2-r )  +  |  BnHn(2,r))e; 


(118) 
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The  total  fields  inside  the  conductor,  that  is,  in  region  II, 
where  a<r<b,  are 


Er  -  iE0eihz  l  in  {  £  [CnJ^(xr)  +  DnH^(xr)]  + 

— oo 

+  I?  fFnJn(xr)  +  GnHn(xr)^ein4>  '  (119) 

E<t,  -  iE»eih2  |  in  {]^(in)[CnJn(xr)  +  V^xn]  + 

-  tFnJA(xr)  +  GnHA(xr>^ein<t>  '  (120) 

E2  -  E0eihz  l  in  |  [CnJn(xr)  +  DnHn(Xr)]ein<t>  ,  (121) 

—  CO 

Hr  =  -iH0eih2  I  in  {  ia  ^I[CnJn(xr)  +  DnHn(Xr)]  + 

■  I  tFnJA(xr)  +  GnHA(xr)^eln<t>  '  (122) 

H4,  “  iH,eihz  l  in  {  £l[CnJ^(xr)  +  DnH^(Xr)]  + 

—  oo 

+  ]^F(in)tFnJn(xr)  +  GnHn(xr)^ein4>'  <123) 

Hz  -  H0eihz  l  in  £  [FnJn(xr)  +  GnHn(xr)]ein<t>.  (124) 


The  total  fields  on  the  interior  of  the  shield,  in  region  III, 
where  r<a,  are  given  by 


Er  =  iE0eihz  J  in  (£  PnJ^(8-r)  +  QnJn(fi<r)  )ein<t>,  (125) 

— oo 

E  «  iE0eihz  l  in  (^(in)  PnJn(fcr)  -  QnJ^lz))ein* ,  (126) 

—  oo 

Ez  «  E0elhz  £  in  £  PnJn(fi-r)ein<1),  (127) 

—  oo 

Hr  =  -iH0elhz  l  in  PnJnUr)  -  |  QnJ^(fi<r))ein<,>,  (128) 

—  oo 

H  “  iH0eihz  l  in  (PnJ^(«-r)  +  j^Un)  QnJn( 2-r )  )ein4>,  (129) 

—  oo 

Hz  =  H„eihz  l  in  £  OnJn(ar)ein4>.  (130) 

—  oo 


To  determine  the  expansion  coefficients  the  boundary  conditions 
must  be  applied  at  the  surfaces  r«a  and  r=b.  The  boundary 
conditions  require  that  the  tangential  components  of  the  E-  and 
the  H-fields  be  continuous  on  the  boundary  surfaces.  The 
condition  that  E^  is  continuous  at  r=b  implies 


-Jn<lb>  +  KIB<1">AnHn'llb>  ‘  BnHA>lb>  ' 

■  E?B,in)[cnJn,xb)  +  DnHn<xb)J  “  tVA(*b)  +  GnHA<*b,J-  <131> 
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The  condition  that  Ez  is  continuous  at  r-b  implies 

fAnV,'b>  '  I  tc„Jn<xb>  +  »nV*b'l- 

The  requirement  that  is  continous  at  r-b  implies 

HE<in>Jn<lb>  +  VA<«»  ♦  - 

£lCnJA<*b>  +  DnHA<*b>]  +  EtE<in'(FnJn^b>  +  GnV*b>]- 

The  requirement  that  Hz  is  continous  at  r-b  implies 

I  V“»  +  I  W*b'  -  £  tpnJn<*b>  +  <W*b>l‘ 

The  requirement  that  is  continous  at  r-a  implies 

sfsUnltCnVxa)  ♦  DnHn(Xa)]  -  [F^lxa)  +  GnH^(xa)  ] 

■  ei,ln>  W*a'  -  VA'l*>' 

The  requirement  that  Ez  is  continous  at  r-a  implies 

I  [Vn(*a>  +  W*a>]  -  |  Vn'*a>- 

The  requirement  that  H.  is  continuous  at  r-a  implies 


(132) 


(133) 


(134) 


(135) 


(136) 
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And,  the  requirement  that  H  is  continuous  at  r-a  implies 

Z 


£  tVn<*a>  +  GnHn<Xa)l 


(138) 


These  are  all  of  the  required  boundary  conditions.  The  full  set 
of  boundary  conditions  constitutes  a  nonhomogeneous  set  of  eight 
equations  in  the  eight  unknowns  An,  Bn,  Cn,  Dn,  F  ,  Gn,  P  ,  and 

Qn<  This  system  may  be  solved  by  the  standard  procedure.  The 

explicit  solutions  for  the  coefficients  is  straightforward  and 
shall  not  be  given  here.  According  to  Kramer's  rule,  the 
solutions  are  given  by  relations  of  the  form 


<|C 

M 

% 

V 

U, 

V  ' 

Cn“ 

u» 

V  ' 

V 

u4 

V  ' 

u, 

V  ' 

V 

u, 

V  ' 

V 

U, 

V  ' 

U. . 

V  ' 

(139) 


where  the  determinants  and  v  are  determined  by  inspection  of 

the  eight  equations  131,  132,  133,  134,  135,  136,  137,  and  138. 
This  completes  the  solutions  for  the  electromagnetic  fields  for 
the  case  of  oblique  incidence.  Due  to  the  fact  that  the  expansion 
coefficients  are  given  by  eight-by-eight  determinants  the 
shielding  effectiveness  for  the  case  of  oblique  incidence  shall 
not  be  calculated. 

Although  the  solutions  of  the  fields  is  somewhat  more 
complicated  than  in  the  case  of  perpendicular  incidence,  the  fact 
that  for  a  good  conductor  only  a  single  mode  is  dominant,  implies 
that  the  behavior  of  the  shielding  effectiveness  should  be  very 
similar  to  the  case  of  perpendicular  incidence.  While  this  is 
almost  certainly  true  for  near  perpendicular  incidence,  for  large 
angles  of  incidence,  and  especially  grazing  incidence,  this  can 
only  be  verified  through  calculation. 


7.  SHIELDING  EFFECTIVENESS  OF  THE  CYLINDER 


The  main  purpose  of  this  investigation  is  to  calculate  and 
plot  the  shielding  effectiveness  verses  frequency  for  hollow 
cylindrical  conducting  shields.  Results  shall  be  presented  only 
for  the  case  of  perpendicular  incidence.  For  perpendicular 
incidence,  the  shielding  effectiveness  is  defined  for  the  TM 
incident  wave,  the  E-polarized  wave,  as  the  ratio  of  the 
magnitude  squared  of  the  E-field  on  the  axis  of  the  cylinder  to 
the  magnitude  squared  of  the  E- field  of  the  incident  wave 


(TM) 


SE, 


"axis 


"inc 


(140) 


Similarly,  the  shielding  effectiveness  for  the  TE  incident  wave 
is  defined  as  the  ratio  of  the  magnitude  squared  of  the  H-field 
on  the  axis  to  the  magnitude  of  the  H-field  of  the  incident  wave 


(TE) 


SE. 


H 


axis 


H 


inc 


(141) 


The  shielding  effectiveness  is  usually  expressed  in  terms  of  dB 
notation.  For  equations  140  and  141  the  shielding  effectiveness 
in  dB  is 


(TM) 


SEe(dB) 


-20  log 


E. 


axis 


"inc 


(142) 


(TE) 


SEh(dB) 


-20  log 


H 

H 


axis 

inc 


(143) 


To  calculate  the  shielding  effectiveness  the  values  of  the  fields 
on  the  axis  must  be  determined.  Using  the  series  expansion  of  the 
E-field  in  region  III  for  the  TM  incident  wave,  equation  29, 
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together  with  the  fact  that  Jn(0)=  1  if  n=0,  and  Jn(0)=  0  if 
n-±l,±2,...,  the  field  on  axis  is  given  by 


(TM) 


Eaxis 


E0d o . 


(144) 


Similarly,  from  the  series  expansion  of  the  H-field  for  the  TE 
case,  equation  47,  the  field  on  axis  is  given  by 

(TE)  Haxis  =  H0D0 .  (145) 


Consequently,  the  shielding  effectiveness  for  the  TM  and  TE 
incident  waves  are 


(TM) 

SEe  =  |de|2. 

(146) 

(TE) 

SEh  «  1  Do | 2 • 

(147) 

The  quantities  d0  and  D0  are  evaluated  in  Appendix  B  and 
Appendix  C,  respectively.  For  the  case  of  perpendicular  incidence 
the  final  results  for  the  shielding  effectiveness  for  the  TM  and 
TE  incident  waves  take  the  following  forms 


(TM) 

se;1/2  - 

l\d  + 

1  Ti*ab' 

(148) 

(TE) 

SE-1/2  - 

|Aa'  + 

1  Tisab' 

(149) 

where 

X  e 

2Ae  * 

sin( 6+i6 ) , 

(150) 

y  * 

2Ae 

cos ( 6+i6 ) , 

(151) 
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a  « 

y*H0(kb)J0(ka)  +  k*H1(kb)J1(ka) , 

(152) 

P  - 

yktJ^kaJH^kb)  -  HQ(kb)  J^ka)  ] , 

(153) 

a'  - 

k*H0(kb)J0(ka)  +  y*H1(kb)J1(ka) , 

(154) 

and 

P’  - 

yk[jQ(ka)H1(kb)  -  J1 (ka)HQ (kb) ] , 

(155) 

where 

A  -  i  |ya|-1/2|yb|-1/2. 

(156) 

and 

6  -  2_1/2Jy J (b-a) . 

(157) 

These  results  were  derived  by  approximating  the  Bessel  and  Hankel 
functions  having  the  complex  arguments  z«=ya  and  z«yb,  by  their 
principle  asymptotic  forms.  Consequently,  these  equations  are  not 
exact,  and  they  do  not  hold  for  all  frequencies.  However,  as 
discussed  in  Appendix  B,  these  equations  are  sufficiently 
accurate  for  all  frequencies  in  the  range  /«oy0a>125 . 

The  shielding  effectiveness  was  calculated  for  a  conducting 
cylinder  having  an  outer  diameter  of  155  mm.,  conductivity 

10  ohm  m  ,  and  thickness  0.1  mm.  The  results  for  the  TM  case 
are  shown  as  a  solid  curve  in  figure  3,  and  the  results  for  the 
TE  case  are  shown  as  the  solid  curve  in  figure  4.  The  low 
frequency  behavior  is  not  very  interesting;  below  10 *  Hz.  the 
shielding  effectiveness  decreases  to  0  dB  almost  as  a  straight 
line  on  a  log-log  plot  like  figure  3.  For  higher  frequencies  the 
shielding  effectiveness  increases  rapidly  as  is  expected  by 
comparison  with  the  behavior  of  a  planar  shield.  Furthermore  a 
fine  structure  appears.  It  is  well  known  that  a  pair  of  parallel 
conducting  planes  gives  rise  to  Fabry-Perot  resonances,  causing 
increased  transmission  of  the  incident  beam,  when  the  distance 
between  the  plates  is  equal  to  integer  multiples  of  a  half 
wavelength.  A  similar  resonant  behavior  is  expected  for  the 
cylindrical  conducting  shell,  causing  an  increase  in  the 
penetration  of  the  cylinder  by  the  incident  wave.  This  is  born 
out  by  the  theory.  As  can  be  seen  in  figures  3  and  4,  the 
shielding  effectiveness  shows  sharp  dips  at  the  resonant 
frequencies  given  by  the  zeros  of  the  equations 
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TM  INCIDENT  WAVE 


0*009  0*006  0*00*  0*006  0*002  0*001  0*0 

(8P)  SS3N3AI103JJ3  0NICH3IHS 
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Shielding  effectiveness  for  the  cylinder:  TM  incident  wave 


Jo(ka)«0,  for  the  TM  incident  wave,  and  (158) 

Jj(ka)-0,  for  the  TE  incident  wave.  (159) 


A  derivation  of  the  resonance  conditions,  equations  158  and  159, 
is  given  in  Appendixes  B  and  C.  The  shielding  effectiveness  drops 
by  about  3  or  4  orders  of  magnitude  at  the  resonance  frequencies 
in  this  example.  Nevertheless,  in  practice,  and  especially  for 
thick  shields  this  effect  is  negligible.  Calculations  of  the 
shielding  effectiveness  for  a  thick  shield  having  a  1.0  mm 
thickness  prove  this.  The  resonances  are  much  less  pronounced  for 
the  thicker  shield  due  to  the  overall  magnitude  of  the 
attenuation  of  the  shield. 

Using  the  solutions  for  the  shielding  effectiveness  developed 
in  Appendix  B  and  Appendix  C  an  approximate  closed  form  result 
for  the  shielding  effectiveness  is  derived  in  Appendix  D.  The 
approximate  formula,  which  is  identical  for  the  TM  and  TE  modes, 
is  given  by 


SE 


(160) 


or,  in  terms  of  dB, 

SEdB  "  10{lo9(  4^-)  -  log(w)  +  (2om0  )1/2(b-a)log(e)  to1/2}.  (161) 


The  shielding  effectiveness  calculated  by  using  equation  161  is 
plotted  as  a  dashed  curve  in  figures  3  and  4.  The  agreement  with 
the  exact  solutions  is  very  good  for  frequencies  in  the  resonance 
regime.  For  lower  frequencies  this  approximation  breaks  down.  As 
is  shown  in  Appendix  D,  this  formula  diverges  as  to  approaches 
zero  and  therefore  cannot  be  used  for  frequencies  below 


«0 


2  1 
y0o  (b-a) * ‘ 


(162) 


However,  the  approximate  formula  161  is  very  useful  since  it 
allows  shielding  calculations  to  be  made  quickly  on  a  hand 
calculator  without  the  need  for  complicated  mathematics  or 
computer  programs. 
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Figure  4.  Shielding  effectiveness  for  the  cylinder:  TE  incident  wave 


8.  ELECTROMAGNETIC  FIELDS  FOR  THE  SPHERE 


A  conducting  sphere  with  inner  radius  r-a  and  outer  radius 
r«b  is  situated  in  empty  space.  The  material  parameters  of  the 
conductor  are  e,  y0,  and  o.  The  conductor  is  assumed  to  be  a  non¬ 
magnetic  material,  so  it  is  possible  to  set  y«y0.  A  linearly 
polarized  planewave  of  frequency  f  is  incident  on  the  sphere.  The 
planewave  propagates  along  the  z-direction  of  a  cartesian 
coordinate  system  with  the  origin  at  the  center  of  the  sphere. 

The  total  field  in  the  region  exterior  to  the  sphere  is  given  by 
the  sum  of  the  incident  field  and  the  scattered  field,  that  is. 


E 


Einc  +  Escat‘ 


The  time  dependence  of  the  field  vectors  is  the  same  as  that  of 

the  incident  wave,  which  is  assumed  to  be  e-ico^,  where  the 

2 

imaginary  unit  i  -  -1  is  used  to  prevent  confusion  with  the 
spherical  Bessel  function  j(z).  The  fields  of  the  incident 
planewave  are 


E 


x 


E.elkz, 

(163) 

H„elkz, 

(164) 

where  k-«/c,  c«l//y0e0  is  the  speed  of  light,  and  H0=/e0/y0  E0 . 
Following  Stratton,  the  fields  may  be  expanded  in  terms  of  the 
vector  spherical  wavefunctions  L,  M,  and  N,  since  these  functions 
form  a  complete,  orthogonal  set  of  solutions  to  the  vector  wave 
eguation  in  spherical  coordinates.  The  procedure  is  well  known. 
The  fields  of  the  incident  planewave  may  be  expressed  by  the 
expansions 


E 


inc 


E,  I  i 

n=l 


n  2n+l 
n(n+l ) 


[M 


o.  In 


iNe, ln^p'kr 


(165) 


H 


inc 


H0e 


ikz 


'nil  n<n+1> 


IK 


e,  m 


1Ni,lnW'<166> 
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where  the  vector  spherical  wavefunctions  are  given  by 


M '  .  ■ 

o,  m 

sine  Vp)PA<cose>  c05*  ®e 

ur 

'  jn(p)  3e~  sin4>  ®<t>' 

(167) 

Ml  1n  * 
e ,  In 

sine  Vp)PA<cose>  sln*  ee 

% 

dP  * 

-  Vp)  35“  cos*  e<J>' 

(168) 

and 

N‘  ,  - 
o ,  in 

jn(P)Pn(cos0)  sin4>  e 

i  dpA 

r  +  p  [pjn(P)]’  gg-  sin<t> 

e0  + 

+  psine  tPVP'J'  pA(cos0) 

cos<J>  e^. 

(169) 

Ni  ,  = 
e.  In 

n  "1  jn(p)P‘(cos8)  cos<J>  e 

i  dpA 

r  +  p  [pjn<P>]'  a?-  cps* 

e0  + 

-  pSIHe  pA,cose)  sln,t'  V 


(170) 


the  prime  denoting  differentiation  with  respect  to  p.  The 
function  jn(P)  is  the  spherical  Bessel  function,  and  the  function 

P^(cosG)  is  the  associated  Legendre  polynomial  with  upper  index 

m*l.  The  wavefunction  L  is  not  required  for  expansions  of  the  E 
and  H  fields  since  the  divergence  of  L  is  not  zero.  The 
wavefunctions  L,  M,  and  N  compose  a  fundamental  set  of  solutions 
of  the  vector  Helmholtz  equation,  in  spherical  coordinates 

-VxV*F  +  k2F  -  0.  (171) 


Furthermore,  the  functions  M  and  N  satisfy  the  relations 

(172) 

(173) 


VxM  -  kN, 
VxN  *  kM. 
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A  derivation  of  the  vector  spherical  wave  functions,  their 
orthoganality  relations,  and  other  properties  is  given  in  [xx]. 
The  scattered  fields  outside  the  sphere  may  be  expressed  in  terms 
of  series  expansions  analogous  to  those  used  for  the  incident 
fields.  For  r>b  the  scattered  fields  are  given  by 


Escat 

Hscat 


«  v  in_2n±l_ 

°n=l  n(n+1) 


V  2n+1 
n(n+l ) 

n=l 


tanMo,ln’  ibnNe,lnV=)cr  ' 
tbnMe,ln+  ianNo, ln^p^kr  ' 


(174) 


(175) 


where  the  wavefunctions  with  the  superscript  3,  M’  and  N* ,  are 
obtained  from  the  defining  relations  for  M1  and  N1  by  replacing 
the  spherical  Bessel  function  jn(P)  with  the  spherical  Hankie 

function  of  the  first  kind  h^(p).  Next,  similar  expansions  for 

the  fields  inside  the  conducting  shell  shall  be  obtained. 

Assuming  a  time  dependence  of  the  form  Maxwell's  equations 

inside  the  conductor  are  given  by 


VxE  ■ 

icoy  0H, 

(176) 

V  *H  = 

(o-icoe  )E, 

(177) 

V'E  = 

0, 

(178) 

V»H  = 

0, 

(179) 

where  the  material  equations  are 


D  = 

eE, 

(180) 

B  - 

y«H, 

dbl) 

J  * 

oE. 

(182) 

The  field  equations  imply  that  the  field  vectors  E  and  H  satisfy 
the  vector  wave  equation 
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0 


(183) 


2 

7xVxE  +  y  E 
VxVxH  +  y2H  *  0, 


where 


iwy0  (o-icoE ) . 


(184) 

(185) 


Therefore,  the  fields  inside  the  conductor,  for  a<r<b,  have  the 
expansions 


E  - 


E-J1inHTH7n-  tcnMo.ln+  Vi, In'  lfnNi,ln-  ' 

(186) 


H  *  -Xh0  y  in-f~T  +  gnM‘  ic  N*  .  +  id  N'  .] 

k  n(n+l )  1  n  e,ln  3n  e,ln  n  o,ln  n  o,inJp»yr 

(187) 


Wavefunctions  of  both  the  first  and  third  kinds  are  required  to 
be  able  to  match  the  boundary  counditions  with  both  the  exterior 
and  interior  fields.  Similarly,  the  fields  inside  the  shielded 
volume,  r<a,  have  the  following  expansions 


E 


H 


p  V  *n  2rvfl 
°nii  n(n+l) 


oo 

-He  l 
n«l 


, n  2n+l 
1  n(n+l) 


^nMo,ln~  iqnNe, ln-lp=kr  ' 
^qnMe,ln+  i^nNo, ln^p*kr  ’ 


(188) 


(189) 


Only  the  wavefunctions  of  the  first  kind  may  be  used  in  these 
expansions  since  the  fields  must  be  finite  at  the  origin. 

The  next  step  in  the  solution  of  the  boundary  value  problem  is 
to  match  boundary  conditions.  The  boundary  conditions  on  the 
field  vectors  require  that  the  tangential  components  of  the  field 
vectors  E  and  H  must  be  continuous  at  the  boundaries  r*a  and  r=b. 
The  condition  that  the  tangential  components  of  E  are  continuous 
at  r*b  Implies 
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(190) 


3n<kb|  +  anhn(ltb)  *  cnhn<',b>  +  dnVvb>- 

J{[pJn(P)]'+  bn[Phn<p|l')p.kb  '  (fn[pVp)J'+  9ntp3n(P)]'}p,Yb- 

(191) 

The  condition  that  the  tangential  components  of  H  are  continuous 
at  r«b  implies 


jn(kb)  +  bnhn(kb)  -  £  [fnhn(yb)  +  (192) 

{[pdnCP)]'  +  antphn(p)l' Wb  "  {cn[Phn(p>]'  +  dJpVp)  1 '  Wb’ 

(193) 

The  condition  that  the  tangential  components  of  E  are  continuous 
at  r«*a  implies 


cnhn<Ya>  +  Vn'**'  '  pnVka>'  (194> 

{fn[phn(p)]'  +  9„[p3n(P)]'}p.ya  -  )r  <3nltP3n<p>]'5p=ka'  <195> 


And  lastly,  the  condition  that  the  tangential  components  of  H  are 
continous  at  r«=a  implies 


J  [fnhn(Ya)  +  gnjn(Ya)]  -  qnjn(ka), 


(196) 


{Cn[phn(p)]'  +  dn[pdn(P)]'} 


p-ya 


Pnttrfn^Jl'^p-ka-  (197> 


These  are  all  the  required  boundary  conditions.  This  set  of 
equations  constitute  a  nonhomogeneous  system  of  eight  equations 
in  eight  unknowns.  The  solutions  of  this  system  yield  the 
expansion  coefficients  for  the  fields.  To  investigate  the 
shielding  properties  of  the  conducting  sphere  it  is  necessary  to 
solve  for  the  fields  inside  the  shielded  volume.  This  implies 
that  it  is  necessary  to  solve  for  the  coefficients  pn  and  q  .  To 

avoid  the  direct  evaluation  of  8*8  determinants  the  system  shall 
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be  reduced  to  a  system  of  4x4  equations  which  shall  then  be 
solved  by  Kramer's  rule.  The  coefficients  a„ ,  b  ,  ,  and 

shall  first  be  eliminated  from  the  eight  equations  to  yield  four 
equations  in  the  four  unknowns  cn,  dn,  ffi,  and  g  .  Having  solved 

the  4x4  system  of  equations,  the  coefficient  pn  may  be  calculated 

from  cn  and  dn  by  using  equation  194,  and  the  coefficient  qn  may 

be  calculated  from  f  and  gn  by  using  equation  196.  To  eliminate 

an,  solve  equation  190  for  an,  and  then  substitute  the  result 

into  equation  193.  Omitting  the  subscript  n  from  all  the 
expansion  coefficients  and  all  the  Bessel  and  Hankel  functions 
for  simplicity,  this  procedure  yields  the  result 

c[rh<Yb>h<kb>  ~  rh(kb)h(yb)]  +  d[ I\  (yb)h(kb)  -  rh(kb)j(yb)]  - 

«  I\ (kb)h(kb)  -  rh(kb)j (kb) ,  (198) 

where  rh(yb)  "  ( [ Phn ( p >  ^  %-yb' 

Fj(kb)  -  {[pjn(P)]'}p,kb/  etc-  <199^ 

The  right-hand-side  of  equation  198  may  be  simplified  further  by 
taking  the  derivative  in  the  r's  yielding 

p[j'(P)h(p)  -  h' (p)j (p) ]p=kb; 

then,  substituting  h(p)  *=  j(p)  +  iy(p),  this  equals 

-ip[3n(p)y A(P)  -  hA<P>V?>Wb- 

The  bracketed  term  is  the  Wronskian  of  jn(P)  and  yn(p),  which 

equals  l/pa.  Thus,  the  left-hand  side  of  equation  198  takes  the 
form 


rj(kb)h(kb)  -  rh(kb)j(kb)  =  j^j. 
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To  recapitulate,  equation  198  takes  the  simplified  form 


c[rh(yb)h(kb)-rh(kb)h(yb)]  +  d^  (yb)h(kb)-rh(kb) j  (yb)  ]  -  jj£. 

(200) 

Next,  to  eliminate  bn  solve  equation  192  for  bn,  and  then 
substitute  the  result  into  equation  191.  This  yields 


f[rh<Yb)h(kb)  -  j£lrh(kb)h(yb)  ]  +  g[l\  (yb)h(kb)  -  ^rh(kb)  j  (yb)  ]  = 

=  Jr  [Tj (kb)h(kb)  -  rh(kb)j(kb)j 

■  $  p[3A'P»hn<pi  -  hA(p,Ve>]p-kb 

■  £  is-  <201> 

To  eliminate  pn,  solve  equation  194  for  pn,  and  substitute  the 
result  into  equation  197.  This  yields  the  equation 


c[rh<Ya>3(ka)  -  (ka)h(ya)  ]  +  d[r_j  (ya)  j  (ka)  -  (ka)  j  (ya)  ]  =  0. 

(202) 

Lastly,  to  eliminate  qn  solve  equation  196  for  qn  and  substitute 
the  result  into  equation  195.  This  yields 

f[rh^ya)d(ka)  "  (ka)h(ya) ]  + 

+  g[lj  (ya) j  (ka)  -  (ka) j  (ya)  ]  -  0.  (203) 
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The  reduced  system  of  equations  is  given  by  the  four  equations 


c[rh(yb)h(kb)  -  rh(kb)h(yb)]  +  d[Tj  (yb)h(kb)  -  rh(kb)j(yb)]  «  ^ 
c[rh(ya) j (ka)  -  (ka)h(ya) ]  +  d[r^ (ya) j (ka)  -  (ka) j (ya) ]  «  0 

f[rh(yb)h(kb)-^*rh<kb)h<Vb)  ]  +  9[rj(Yb)h(kb)-jJrh(kb)j(yb)]« 

f[rh<Ya^  (ka)-j£ar_j  (ka)h(ya)  ]  +  g^  (ya) j  (ka)-j^r_j  (ka)j  (ya)  ]  -  0. 


The  matrix  of  coefficients  has  the  convenient  block  form 


D 


W, 

0  0 
0  0 


0  0 
0  0 

W,  ' 


(204) 


where  wx  = 


[rh(yb)h(kb)  -  rh(kb)h(yb)3  [r^ (yb)h(kb)  -  rh(kb)j(yb)] 
[rh(ya)j(ka)  -  r_j  (ka)h(ya)  ]  (ya)  j  (ka)  -  (ka)  j  (ya)  ] 


and  W,  *= 

[rh(yb)h(kb)  -  ^rh(kb)h(yb)]  (yb)h(kb)  -  ^*rh(kb)  j  ( yb)  ] 

[rh(ya)J(ka)  -  (ka)h(ya)  ]  [r^  (ya)  j  (ka)  -  (ka)  j  (ya )  ] 


The  determinant  D  is  given  by  D  -  W,wt.  Evaluating  the 
determinants  w,  and  W,,  and  rearranging  terms  to  group  together 
the  functions  containing  the  complex  arguements  z-ya  and  z*yb 
gives  the  result 


and 


W,  -  j<ka)h(kb)[rj(ya)rh(yb) 
+  (ka)h(kb)[h(ya)r.j  (yb) 
+  j (ka)Th(kb)[rh(ya)j (yb) 
+  Tj (ka)Th(kb) [j (ya)h(yb) 


rj(Yb)rh(Ya>^l  + 
j (ya)rh(yb)]2 

rj (ya)h(yb)]3 

h(ya)j (yb) ]4, 


w2  =  j(ka)h(kb)[rj(ya)rh(yb)  -  rh(ya)T. (yb) ]x  + 
+  (ka )h(kb)j^a  [h(ya )  ( yb)  -  j  (ya)rh(yb)  ]2 

+  j(ka)rh(kb)^I[rh(ya)j(yb)  -  (ya)h(yb) ]3 
+  (ka)rh(kb)|*[j (ya)h(yb)  -  h(ya) j (yb) ]4 . 


(205) 


(206) 


The  bracketed  terms  have  been  numbered  from  1  to  4,  for  further 
reference.  Now  it  is  possible  to  solve  for  the  coefficients  c  and 
d  to  obtain  p.  Using  Kramer's  rule,  the  coefficient  c  is  given  by 
c  -  u i /D ,  where  » 


Ui 


j~5  [iyybjhfkb)  -  rh(kb)j(yb)] 
0  [ Tj (ya) j (ka )  -  (ka) j (ya) ] 

£  EE  0 

0  0 


0  0 
0  0 


W, 
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thus , 


E5  [r^ (ya)j (ka)  -  r^(ka)j(ya)]  W4. 


(207) 


Similarly,  solving  for  d,  d  =  Ua/D,  where 

[rh(yb)h(kb)  -  rh(kb)h(yb)]  ^ 

[ r_.  (ya)j  (ka)  -  (ka) j  (ya)  ]  0 

Ua 

t 

0  0 


-i 

JcE 


0  0 
0  0 

wa 


thus 


£5  [rh(ya)j(ka)  -  ^ (ka)h(ya) ]  wa 


(208) 


Now  solve  for  p  from  equation  194.  The  solution  for  p  is 
P  -  p[h'(P)j(P)  -  V  (PlhlPlJp^aW,'1  -  (209) 


Next,  solve  for  f  and  g  to  obtain  q.  The  solution  for  f  is  given 
by  f  -  U,/D,  where 


U, 


W, 

0  0 

0  0  [rj(yb)h(kb)  -  ^Irh(kb)  j  (yb)  ] 

0  0  0  [r^  (ya)j  (ka)  -  (ka) j  (ya)  ] 
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thus  f  - 


£  EE  [Tj  <ya) J  (ka)  -  (ka)  j  (ya)  ]  Wt .  (210) 


The  solution  for  g  is  g  «  U4/D,  where 


U4 


W, 


0 

0 


-i 

EE 


0  0  [lj(yb)h(kb)  -  ]^Irh(kb)  j  (yb)  ]  £ 

0  0  [r^  (ya)j  (ka)  -  £1^  (ka)  j  (ya)  ]  0 


thus  g  -  £  Jg  [rh(ya)j(ka)  -  £1^ (ka)h(ya) ]  w,.  (211) 


Now  solve  for  q  from  equation  196.  The  solution  for  q  is 


q 


y*  +i  +i 
k*  EE  ya 


(212) 


The  solutions  fu*.  the  expansion  coefficients  which  have  been 
obtained  in  this  section,  and  p  and  q  in  particular,  are  exact. 
Although  the  integer  subscript  n  has  been  omitted  (p=pn,  q*=qn, 

etc.),  the  results  hold  for  all  n»l,2,3,....  Thus,  the  exact 
solution  for  the  electromagnetic  fields  for  scattering  from  a 
spherical  conducting  shell  have  been  obtained. 

Next,  the  explicit  solutions  for  p  and  q  shall  be  used  to 
obtain  the  shielding  effectiveness  of  the  sphere.  Since  the 
shielding  effectiveness  only  depends  on  the  cavity  fields,  the 
fields  in  the  conducting  walls  and  in  the  exterior  region  are  not 
required.  Consequently ,  explicit  expressions  for  the  coefficients 
a,  ’  ,  c,  d,  f  and  g  are  not  required  in  the  analysis. 
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9.  SHIELDING  EFFECTIVENESS  OF  THE  SPHERE 


The  shielding  effectiveness  shall  be  calculated  with 
reference  to  the  fields  at  the  center  of  the  sphere.  The 
shielding  effectiveness  is  most  commonly  defined  as  the  ratio  of 
the  magnitude  squared  of  the  electric  field  at  the  center  of  the 
sphere  to  the  electric  field  of  the  incident  planewave. 


SE 


e 


E( r«0 ) 


'inc 


(213) 


This  definition  shall  be  called  the  electric  shielding 
effectiveness.  Similarly,  the  shielding  effectiveness  may  be 
defined  as  the  magnitude  squared  of  the  ratio  of  the  magnetic 
field  at  the  center  of  the  sphere  to  the  magnetic  field  of  the 
incident  wave. 


H( r=0 ) 


H 


inc 


(214) 


This  definition  shall  be  called  the  magnetic  shielding 
effectiveness.  Another  definition  may  be  given  in  terms  of  the 
energy  of  the  fields  as  the  ratio  of  the  energy  density  at  the 
center  of  the  sphere  to  the  energy  density  of  the  incident  wave 


SE 


?  U ( r=0 ) 

u4 

inc 


(215) 


This  definition  shall  be  called,  simply,  the  electromagnetic 
shielding  effectiveness.  The  factor  of  two  is  included  so  the 
shielding  effectiveness  goes  to  1  as  the  frequency  f  goes  to 
zero.  It  will  be  shown  that  for  planewave  incidence 


SE  =  SEQ  +  SE  . 

e  m 


(216) 
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In  terms  of  the  complex  fields,  and  in  free  space,  the  energy 
density  is  given  by  the  relation 


U  -  ^[e0E»E*+ 

(217) 

The  energy  density  of  the  incident  planewave  is 

12  12 
u„  -  ±e0e;  - 

(218) 

To  proceed  with 
first  note  that 

the  calculation  of  the  shielding  effectiveness, 

since  1„(r«0)  -  0  for  n-1,2,3...,  then 

Jn 

^Mo,m^p-o  *  °'  n-i,2,3.... 

(219) 

(Me,ln)p«0  “  °'  n«l,2,3. . .  . 

(220) 

Also,  from  the  identity 


[P3n(p)]'  *  2^1^  (n+1)jn-l(p)  "  nJn+l (p>  ] '  n-i  ,2,3 - 


and  the  fact  that 


Lim  p^ 
p->0  p 


1/3  for-n-1, 

0  for  n-2,3,4. . . , 


the  only  nonzero  wavefunction  at  p-0  of  the  type  N  is 


(N*  nlp.Q  *  §  [sin0sln<t>  ®r  “  cosesin*  eQ  +  cos*  e^],  (221) 


(Ng  . .  ] p_Q  -  -3  [cosecos*  e0  +  sin*  e^]. 


(222) 
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Therefore,  the  fields  take  on  a  particularly  simple  form  at  the 
center  of  the  sphere  since  only  one  term  in  the  series  expansions 
is  nonzero,  the  fields  are  given  by 

E  '  2  E°qiNe,ll lp=0' 

H  =  2  H°^lNo, 11 1 p=0 ' 

Therefore,  the  magnitudes  of  the  electromagnetic  fields  at  the 
origin  are 


(223) 

(224) 


E  |  2  - 

1 e0 1 2 1 q i 

I2 

(l)2 

1  p*0 

lEo 1 2 1 q i I2, 

(225) 

H  |  2  = 

|h0|2Ip. 

I2 

If)2 

No<ll"No,ll 

lp=0 

lHo 1 2 1 p i I2. 

(226) 

The  shielding  effectiveness,  then,  is  given  by 


SE  = 
e 

|g. I2. 

(227) 

SE  = 
m 

Ip.!2. 

(228) 

se  -  |Pl 

I2+  k.l2. 

(229) 

This  establishes  the  relation  between  the  three  definitions  of 
shielding  effectiveness  given  by  equation  216.  Substituting  for  p 
and  q  from  equations  209  and  212,  the  final  formulas  for  the 
electric  and  magnetic  shielding  effectiveness  take  the  form 


SEe1/2  *  t|t  "  i*/y|(ka)(kb)|w2|, 

SEm1/2  =  JpT  “  jya|(kb)|w,|. 


(230) 

(231) 
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where  W2  and  W,  are  given  by  equations  205  and  206,  respectively. 
Similarly,  the  electromagnetic  shielding  effectiveness  takes  the 
form 

SE  -  ( |k/y | (ka) (kb) |w, | )~2+  ( |ya | (kb) |wx | )~2  (232) 


In  terms  of  dB,  the  shielding  effectiveness  takes  the  form 


SEe  (dB)  -  -20  LoglqJ  -  20  Log(  |k/y  |  (ka)  (kb)  |wa  |  ) ,  (233) 


SEm  (dB)  -  -20  Log JPl |  -  20  Log ( | ya | (kb) | w, | ) ,  (234) 

and 

SE  (dB)  =  -10  Log( |Pl |2+  |q, |2) 


«  -10  Log{( |k/y|(ka)(kb) |Wa|)~2+  ( |ya | (kb) |w, | )'2} .  (235) 


All  these  results  are  exact.  The  spherical  Bessel  and  Hankel 
functions  which  are  needed  to  calculate  these  quantities  are 
given  in  Appendix  F.  Appendix  F  also  gives  various  approximations 
for  the  spherical  Bessel  functions  which  are  used  in  the 
numerical  calculations.  For  the  low-frequency  calculations  the 
quantities  pj  and  q,  are  calculated  from  equations  209  and  212  by 
calculating  w,  and  W,  directly  from  equations  205  and  206.  For 
the  high-frequency  calculations  the  brackets  [ j ^ ,  []2,  []3,  and 

[]4  which  appear  inequations  205  and  206  cannot  be  evaluated  term 

by  term  for  reasons  explained  in  Appendix  F.  Therefore  some 
approximations  must  be  made  to  evaluate  these  quantities.  The 
details  of  these  approximations  are  given  in  Appendix  F. 
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10.  NUMERICAL  CALCULATIONS  FOR  THE  SPHERE 


The  shielding  effectiveness  is  now  calculated  for  a  spherical 
shell  with  inner  radius  a*77.4mm  and  outer  radius  b=77.5mm.  The 

7  —1—1 

conductivity  of  the  sphere  is  o-io  ohm  m  which  is 
representative  of  good  conducting  materials  such  as  copper  or 
aluminum.  The  results  of  the  numerical  calculations  for  the 
electromagnetic  shielding  effectiveness  at  low  frequencies  is 

3 

shown  in  Figure  5.  This  plot  covers  the  frequency  range  from  10 

to  10  Hz.  As  the  frequency  approaches  zero  the  shielding 
effectiveness  goes  to  one,  that  is,  zero  dB,  as  expected.  As  the 
frequency  goes  to  zero  the  magnetic  shielding  effectiveness  SEm 

approaches  one  since  a  static  magnetic  field  is  not  shielded  by 
the  conductor.  On  the  other  hand  the  electric  shielding 
effectiveness  SEe  goes  to  zero  since  a  static  electric  field 

cannot  penetrate  the  conductor,  that  is,  the  electric  field 
inside  the  cavity  is  zero.  Therefore,  the  electromagnetic 
shielding  effectiveness  goes  to  l  since  SE  -  SEe+  SEm. 

3 

It  is  remarkable  that  for  frequencies  from  about  7*10  to 

7 

1*10  Hz  the  shielding  effectiveness  is  practically  a  straight 
line.  The  slope  of  this  part  of  the  graph  is  20  dB  per  octave.  A 
straight  line  approximation  is  roughly 


SE  (dB)  -  201og(f)  +  C, 


where  the  point  at  which  the  straight  line  intersects  SE=0  dB  is 

3 

C  -  -201og(f0),  where  f0=  5*10  Hz.  Therefore  the  shielding 
effectiveness  is  approximately  given  by 


SE  (dB)  «  201og(f/f o ) , 


(236) 


or 


SE  -  (f/f0)"2. 


(237  ) 


This  is  a  very  interesting  result.  According  to  equation  237  the 
shielding  effectiveness  goes  as  one  over  the  frequency  squared. 
Since  the  shielding  effectiveness  is  equal  to  the  square  of  the 
transfer  function  of  the  cavity  fields  with  respect  to  the 
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ELECTROMAGNETIC  SHIELDING  EFFECTIVENESS 
a-77A  mm. 
b=77.5  mm. 


Shielding  effectiveness  for  the  sphere:  low  frequency  results 


incident  fields  this  simple  result  suggests  that  analytical 
expressions  may  be  obtainable  in  the  case  of  an  incident 
electromagnetic  pulse.  This  idea  will  be  developed  momentarily. 
First  the  low  frequency  results  for  the  electric  field  shielding 
effectiveness  will  be  discussed.  No  graphical  results  are  given 
for  this  case  since  they  are  not  very  interesting.  The  numerical 
results  show  that  the  electric  shielding  effectiveness  decreases 
approximately  as  a  straight  line  over  at  least  seven  orders  of 

7 

magnitude  from  1  Hz  to  10  Hz.  The  slope  of  the  straight  line  is 
-20  dB  per  octave.  Therefore,  the  electric  field  shielding 
effectiveness  is  given  by  the  approximate  equation 


SEe(dB) 


-201og(f)  +  C. 


(238) 


At  f=l  Hz,  the  calculations  give  SEg=C=  283.8  dB,  and  therefore 


SEe(dB)  =  -20log( f/f , ) , 


(239  ) 


se6  -  (f/f,r. 


(240) 


where  f, =1.549x10  . 

The  approximate  forms  for  the  shielding  effectiveness  just 
obtained  imply  the  relations 


m°± 

Hinc 


(w/co o  )_1, 


(241) 


HOI 

E  •  _ 


(w/w, ) +1 . 


(242) 


These  simple  transfer  functions  may  be  used  to  obtain  an 
analytical  expression  for  the  fields  at  the  origin  when  the 
incident  planewave  is  a  pulse.  In  cases  when  the  frequency 
spectrum  of  the  pulse  is  bandlimited  to  the  range  of  validity  of 
equations  241  and  242  the  response  may  be  obtained  by  using 
Laplace  transforms.  For  example,  in  the  case  of  a  high  altitude 
nuclear  burst  the  resulting  electromagnetic  pulse  will  have  a 
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waveform  which  is  approximated  very  well  by  the  double 
exponential 


Ei(t)  -  E„y  (e"pt-  e~at],  t>0. 


(243) 


8—1  6—1 
where  a«4. 76*10  s  is  the  rise-constant,  p=4*10  s  is  the 

decay-constant,  Eo-50  kv/m  is  the  peak  electric  field,  and  y=1.05 

is  a  normalization  factor.  For  the  sphere  diameter  chosen  in  this 

report  the  frequency  spectrum  of  this  pulse  lies  almost  entirely 

within  the  range  of  validity  of  241  and  242.  Therefore,  using  241 

and  242,  the  attenuated  fields  at  the  center  of  the  sphere  will 

have  the  Laplace  transforms 

H(S)  =  |  H±(S),  (244) 

E(S)  =  BS  E^S),  (245) 


where  E^s)  and  H^s)  are  the  Laplace  transforms  of  the  incident 

pulse,  A=2nf 0-3.14*104,  and  B=l/w, «1 . 03*10~15 .  From  243  the 
Laplace  transform  of  the  incident  electromagnetic  pulse  is 


Fi(s>  ■  (si^pr  <246 

where  F  stands  for  either  H  or  E.  Therefore,  taking  the  inverse 
transform  the  fields  at  the  center  of  the  sphere  are  found  to  be 


H(t) 


AH,y  [£(l  +  e~at) 


+  e 


-Pt 


)]. 


(247  ) 


E ( t )  *  BE0y  [ae~at-  pe-pt]. 


(248) 


The  analytical  solution  for  E(t)  is  not  of  much  practical  value 
because  the  attenuation  factor  A  is  so  great.  Nonetheless,  these 
expressions  are  a  good  approximation  for  the  EM P  response  of  a 
spherical  shield  of  these  dimensions,  i.e.,  a=77.5mm.  As  a  final 
comment,  note  that  it  is  not  necessary  to  use  Laplace  transforms 
directly.  Using  the  properties  of  the  Laplace  transform  the 
expressions  244  and  245  imply  that  the  response  in  the  time 
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domain  H(t)  and  E(t)  are  just  the  integral  and  the  derivative  of 
the  respective  inputs  H^(t)  and  E^(t). 

Consider  next  the  high-frequency  results  for  the  spherical 
shield  shown  in  Figure  6.  These  high  frequency  calculations  were 
made  by  making  certain  approximations  to  the  spherical  Bessel  and 
Hankel  functions  which  are  detailed  in  Appendix  F.  The  results 
which  are  given  in  Figure  6,  however,  are  entirely  accurate.  The 
most  prominant  feature  in  this  Figure  are  the  sharp  dips  above 

Q 

10  Hz.  These  are  resonances  similar  to  those  obtained  for  the 
cylindrical  shield.  The  resonance  condition  for  the  spherical 
shield  are  derived  in  Appendix  F.  The  resonances  occur  when 
either  of  the  following  conditions  is  met: 


0, 

(249) 

0. 

(250) 

The  first  few  zeros  of  equation  250  are  listed  in  Appendix  F.  In 
general,  the  shielding  effectiveness  is  reduced  by  three  or  four 
orders  of  magnitude  at  resonance.  However  this  effect  is 
negligible  when  compared  to  the  overall  magnitude  of  the 
attenuation  at  these  frequencies.  As  in  the  case  of  the 
cylindrical  shield  an  approximate  formula  for  the  shielding 
effectiveness  is  derived  in  Appendix  F  which  has  the  form 


SE 


8ooe  q  e~/2oooy  c  (b-a ) 


(251) 


This  approximate  formula  is  plotted  as  a  dashed  curve  in 
Figure  6.  Note  that  this  is  the  same  formula  obtained  for  the 
cylinder,  equation  161,  except  for  a  factor  of  2.  A  comparison 
between  this  formula  and  the  exact  solution  shows  that  this  is  a 
very  good  approximation  for  frequencies  in  the  resonance  regime. 
As  explained  in  Appendix  F,  this  formula  has  a  lower  limit  on  its 
range  of  application  since  it  diverges  for  lower  frequencies. 
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Shielding  effectiveness  for  the  sphere:  high  frequency  results 


CONCLUSIONS 


The  theory  of  electromagnetic  shielding  has  been  presented  for 
a  circular  conducting  cylinder  and  a  conducting  sphere.  Starting 
from  first  principles,  the  complete  solutions  for  the 
electromagnetic  fields  were  developed  first  for  the  scattering  of 
a  planewave  by  a  cylindrical  conducting  shell,  and  second  for  a 
spherical  conducting  shell.  The  shielding  effectiveness  of  the 
cylinder  was  defined  with  respect  to  the  fields  at  the  center  of 
the  cylinder  (on  axis),  and  the  shielding  effectiveness  of  the 
sphere  was  defined  with  respect  to  the  fields  at  the  center  of 
the  sphere.  In  each  case  the  fields  at  the  origin  have 
a  particularly  simple  form  since  all  the  terms  in  the  series 
expansion  vanish  except  for  the  first  one,  and  therefore  it  is 
possible  to  exibit  a  single  expression  for  the  shielding 
effectiveness.  This  is  desirable  from  a  practical  point  of  view 
since  problems  with  convergence  of  infinite  series  may  be  avoided 
in  the  numerical  calculations. 

The  most  outstanding  feature  of  the  shielding  effectiveness  is 
the  existence  of  resonances.  In  general,  a  series  of  discrete 
resonances  occur  for  high  frequencies  where  ka>l.  The 
electromagnetic  resonances  for  the  cylinder  and  the  sphere  have 
been  investigated  in  detail  and  the  resonance  conditions  have 
been  derived  for  both  structures.  In  addition,  an  asymptotic 
formula  for  the  shielding  effectiveness  has  been  developed  which 
appears  to  be  a  good  approximation  in  the  resonance  regime.  To 
this  authors  knowledge  this  result  is  new  and  it  has  not  been 
published  in  the  open  literature.  While  this  formula  does  not 
reproduce  the  resonances,  it  does  account  for  the  average 
magnitude  of  the  shielding  effectiveness  in  the  high  frequency 
limit. 

Numerical  calculations  were  performed  for  both  the  cylinder 
and  the  sphere,  and  the  shielding  effectiveness  was  plotted  as  a 
function  of  frequency  for  a  shield  with  an  outer  diameter  of 
155  mm,  a  thickness  of  0.1  mm,  and  a  conductivity  of 

107  ohm  1m  1.  The  results  are  given  in  figures  3,  4,  5,  and  6.  As 
expected,  the  features  of  the  shielding  effectiveness  are 
qualitatively  and  quantitatively  similar  for  both  the  cylinder 
and  the  sphere.  The  magnitude  of  the  frequency  dependence  of  the 
shielding  effectiveness  is  the  same  for  both  structures.  Both 
structures  show  similar  resonant  behavior,  although  not  at  the 
same  resonant  frequencies.  And,  both  structures  have  identical 
asymptotic  forms  for  frequencies  in  the  resonance  regime.  The 
almost  linear  increase  of  the  shielding  effectiveness  at  low 
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frequencies  which  was  discovered  for  the  sphere  could  not  be 
demonstrated  analytically.  However,  this  is  believed  to  be  a 
general  feature  of  the  solutions  for  both  the  cylinder  and  the 
sphere . 


APPENDIX  A 

NUMERICAL  CALCULATION  OF  BESSEL  FUNCTIONS 
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To  calculate  the  series  expansion  coefficients  it  is  necessary 
to  evaluate  Bessel  and  Hankel  functions  of  both  real  and  complex 
arguments.  Specifically,  it  is  necessary  to  evaluate  J0(z), 

J i ( z ) ,  H0(z),  and  Hj(z),  and  the  first  derivatives  of  these 
functions,  where  z  is  in  general  complex.  The  Bessel  function  of 
the  first  kind  is  defined  by  the  series 


Jn<2> 


oo  k 

y  (r1) 

k"0  k! (k+n) ! 


(f)2k+n 


n*0 ,1,2,... 


(A-l) 


The  Neumann  function,  or  Bessel  function  of  the  second  kind,  is 
defined  here  for  positive  integer  indicies  by 


vz> 


n[Y  +  lJn(z) 


1 

Tt 


n-1 

X 

k*=0 


(n-k-1 ) ! 
— Jc!" 


(|)2k'n 


+ 


1  “  (-l)k  r  z  *2k+n 

i  klQ  kl (k+n) !  I2J 


{l+i+...+£  +  1+1+...+^} 


(A-2) 


In  the  case  n=0,  the  second  term  in  this  definition,  the  finite 
sum  from  k=0  to  n-l  is  omitted.  This  term  is  not  present  in  the 
definition  of  Y0(z).  For  negative  indicies  the  Bessel  functions 
shall  be  defined  by 


J-n<z> 

=  (-l)n  Jn(z)  , 

(A-3) 

y-n<z> 

=  (-l)n  Yn(z)  . 

(A-4) 

The 

by 

Hankel  functions  are  defined  in  terms  of  the  Bessel 
the  relations 

functions 

HA(Z)  * 

Jn(z)  +  iYn(z) 

(A-5) 

H*(z)  = 

Jn(z)  "  iYn(z) 

( A-6  ) 
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The  Bessel  functions  are  calculated  in  two  different  ways 
depending  on  the  magnitude  of  z.  The  series  expansions  A1  and  A2 
are  used  when  |z|£  10,  and  the  asymptotic  series  expansions  are 
used  when  |z|>10.  Using  the  defining  equation  A1  the  real  and 
imaginary  parts  of  the  Bessel  functions  are 


Re[JnU)]  -  Jo  (f)2lt+n  cos[(2k+n)6)  ,  (A-7) 


00  K  * 

Im[Jn(z)]  =  l  k» (k^A)!  (f)2k+n  sin[(2k+n)0]  ,  (A-8) 

k=0 


i8 

where  z  =  re  .  These  series  are  used  to  calculate  the  Bessel 
functions  J„(z)  and  J,(z)  when  |z|£l0.  If  the  series  is  truncated 
after  the  26th  term,  then  the  error  is  of  the  order  of  the  first 

term  neglected,  roughly  10-16  for  |z|=10.  For  the  Bessel  function 
of  the  second  kind  the  real  and  imaginary  parts  are  given  by 


Re [ Y 0 ( z ) ]  -  f  {Re[j0 (z) ] (y+ln{r/2) )  -  eim[j0(z)]  -  S0)  , 

(A-9) 

Im[ Y0 ( z ) ]  =  |  {lm[ J0 (z ) ] ( Y+ln( r/2 ) )  +  6Re[j0(z)]  -  T0}  , 

(A-10) 

He[Y,(z)]  =  |  [Re[  J  j  ( z  )  ]  (y+ln(  r/2 )  )  -  eim^tz)]  -  -pcose  -  S,}  , 

(A- 11 ) 

Im[  Y,  ( z )  ]  -  ~  {lm[  J ,  ( z )  ]  ( y+ln  ( r/2 ) )  +  SRefj^z)]  +  |sine  -  T.)  , 

( A-12 ) 


wnere  y  -  0.577215664902  is  Euler's  constant;  and  S„,  S,,  T0, 
and  T,  are  given  by  the  following  series: 
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( A— 1 3 ) 


m  k  2k 

l  ^7^  r  ?  (l+|-+i+. .  .+7)  cos ( 2k6 )  , 

-1  22K(k! r  2  J  K 


To  «  J  (l+i+i+.  •  .+7)  Sin(2k6)  , 

k«l  22  (k ! )l  2  i  K 


(A-14) 


•  ](  2k+l 

S,  -  E  4t'M  — -  [2(l+4+i+---+7)  +  ttt]  COS(2k+l)6,  (A-15) 

k*0  22K+2k!  (k+1) !  2  J  *  K+1 


k  _2k+l 


T,  -  l  -C-U  -I -  [2(l+i+i+. .  .+7)  +7^-]  sin( 2k+l ) e .  (A-16) 

k=0  22K+2k! (k+1)  •  2  J  k  K+i 


To  calculate  Y0(z)  and  Yj(z)  these  series  are  truncated  after  the 
26th  term. 

If  the  magnitude  of  z  is  greater  than  10,  |z|>10,  then  the 
Bessel  functions  are  calculated  by  using  the  asymptotic  expansion 


Jn<2> 


V2> 


where 


V2> 


(f)1/2  z"1/2{Pn(z)cos[z-(2n+l)J]  -  0n(z)sin[z-(2n+l)J]}, 

( A-l 7 ) 

(|)1/2  2-1/2{Pn(z)sin[z-(2n+l)J]  +  Qn( z )cos[ z- ( 2n+l  )J] } ; 

\ 

( A-l 8 ) 


,  ( 4n2-l ) ( 4n2- 9 )  .  (4n2-l)(4n2-9)(4n2-25)(4n2-49) 

1  -  5  +  4  ”  *  *  • ' 

2!(8z)  4!(8z) 


0n(z) 


( 4n2-l )  (4n2-l) (4n2-9) (4n2-25)  , 

82  "  3!(8z)3 


(A-19 ) 
(A-20) 
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The  cooresponding  asymptotic  expansion  for  the  Hankel  function  of 
the  first  kind  is 


HA  <z> 


-  (|)l/2  2'1/2  {Pn(®>  +  iQn ( z ) )  ( A-21 ) 


When  4  terms  are  retained  in  equations  A-19  and  A-20  the  error  in 
the  calculated  values  of  the  Bessel  functions  j0(z),  J,(z), 

Yo(z),  Y,(z),  or  the  Hankel  functions  H0(z),  or  H,(z),  is  less 

than  1*10  for  |z|>10.  The  accuracy  improves  as  |z|  increases. 

The  derivatives  of  the  Bessel  functions  are  calculated  from 
the  functions  themselves  by  means  of  the  following  identities: 


af1  --*.<*>  -  ; 


( A-22 ) 


JA(2) 


Jn-1(2)  2  Jn{2>  '  n-1,2,3...  ; 


( A-23 ) 


yA(2) 


Yn-l^2^  2  Yn^z^  '  0*1,2, 3...  ; 


(A-24) 


HA(z) 


Hn_i ( 2 )  "  §  Hn(z)  ,  n-1,2,3 _  ; 


( A-25 ) 


where  Hn(z)  stands  for  either  H^(z)  orH^(z),  the  Hankel  function 
of  the  first  or  second  kinds. 
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APPENDIX  B 
CALCULATION  OF  d 0 
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For  the  physical  parameters  considered  in  this  study,  and  for 
the  frequency  range  of  interest,  the  Bessel  functions  may  be 
evaluated  for  the  complex  arguments  z«ya  and  z=yb,  by  using  the 
principle  asymptotic  forms  for  the  Bessel  and  Hankel  functions. 
The  asymptotic  form  for  the  Bessel  function  is  given  by 

Jn(z)  =  (f)1/2  z~1/2  cos[z-(n+|)f]  ' 

=  (-^—]1/2  {eitx~^n+2^2^  e_y  +  e-i[x-(n+^)^]  e+y]  ;  (B-l) 


and,  for  the  Hankel  function  of  the  first  kind  by 
Hn(z)  =  (f)1/2  z"1/2 

=  (|)1/2  z~1/2  eifx_^n+2^2^  e~y  ;  (B-2 

where  z  «  x+iy.  A  further  approximation  is  helpful.  Using 
equation  21,  and  considering  the  relevant  physical  parameters, 

that  is,  f  £  5*10 4  Hz,  o=107  ohm_1m_1,  and  b>a  £  7.5  mm,  yields 

| Y |  =  ( wy Oo) 1/2  >  1985, 

and  Arg(y)  =  \  tan-1^)  .  J  , 

which  together  imply 

y  =  2“1/2|ya|  >  10. 

Therefore,  the  first  term  in  equation  B-l  may  be  neglected 
compared  to  the  second  term.  This  yields  the  further 
approximation 


The  calculation  of  d0  involves  the  calculation  of  a  numerator 
and  a  denominator.  The  denominator  is  given  by  equation  41,  with 
n*0.  Evaluating  the  determinant  yields 

V  -  yH0(kb)  [  yJQ(ka)  ( J’ (yb)H£ (ya )  -  J£(ya)H£(yb) )x  + 

kJ'(ka)  (J0(ya)H‘ (yb)  -  J' (yb)H0(ya) )2  ]  + 

-kH'(kb)  [  kJ^(ka)  ( JQ (ya)HQ (yb)  -  JQ (yb)HQ (ya) ) 3  + 

+  yJQ(ka)  ( JQ (yb)H^ (ya )  -  J£(ya)H0(yb) )4  ]  .  (B-4) 

The  inner  brackets  have  been  numbered  from  1  to  4.  These  inner 
brackets  may  be  evaluated  by  first  using  the  identities  A-21  for 
the  derivatives,  and  then  substituting  the  required  asymptotic 
forms  B-2  and  B-3.  The  evaluation  of  bracket  number  1  according 
to  this  procedure  goes  as  follows: 

(  )1  *  J£<yb)H’(ya)  -  J’(ya)H£(yb) 

-  J1(yb)H1(ya)  -  J^yajH^yb) 

=  AC[e~i(Xa'x‘  >  c(y*"y,)  -  ei(x,-xa)  e-<y»-yi)]f 

where 


II 

< 

|ya|-1/2  |yb| 

-1/2 

9 

c  -  V  \ 

( B— 7  ) 

ya  - 

X| +iy , ,  and 

yb  - 

x.+iy, . 

Using  the  approximation:  Arg(y)  «  tt/4,  which  is  valid  for  all 
the  physical  parameters  of  interest,  implies 

x,-x,  -  y,-y,  *  6,  where  6  **  2"  |y|(b-a).  (B-8) 


(B-5) 


(B-6  ) 
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Using  this  result,  equation  B-8,  the  first  bracket,  equation  B-6, 
takes  the  following  form 

(  )x  «  AC{e"i6e6-  ei6e"5J}  (B-9) 

=  A/2  {  [cos ( 6 )sinh{ 6 )-sin( 6 )cosh( 6 ) ]  + 

-  i[cos ( 6 ) sinh( 6 )+sin( 6 )cosh( 6 ) ]  )  .  (B-10) 

This  latter  form  is  useful  for  making  calculations.  Evaluating 
the  next  bracket,  bracket  number  2,  in  exactly  the  same  way 
proceeds  as  follows: 

(  )2  =  J0(ya)H£(yb)  -  J£(yb)H0(ya) 

J1(yb)H0(ya)  -  JQ ( ya)^ (yb) 

-i6  6  i6  — 6 i 
-  ACile  e  +  e  e  } 

=  a/2  {  [cos ( 6 )cosh( 6 )+sin( 6 )sinh( 6 ) ]  + 

+  i[cos{6 )cosh(6 )-sin(6 )sinh(6 ) ]  }.  (B-13) 

This  completes  the  evaluation  of  bracket  number  2.  Proceeding 
similarly  with  bracket  number  3: 

(  )3  =  J0(ya)H0(yb)  -  H0(ya)J0(yb)  (B-14) 

=  AC{ei6e"6  -  e'l6e6}  (B-15) 

*  a/2  {  [ sin( 6 )cosh( 6 )-cos ( 6 )sinh( 6 ) ]  + 

+  i[cos ( 6 ) sinh( 6 ) +sin( S )cosh( 6 ) ]  ).  (B-16) 


( B— 11) 


(B-12 ) 
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This  completes  the  evaluation  of  bracket  number  3.  Lastly,  the 
fevaluation  of  bracket  number  4  goes  as  follows: 

(  )4  -  J0(yb)H'(ya)  -  J£(ya)H0(yb)  (B-17) 

-  J^yajH^yb)  -  H^yaJJ^yb) 

-  ACi{ei6e-6  +  e“i6e6)  (B-18) 

*  a/2  {  [cos(6)cosh(6)+sin(6)sinh)6) ]  + 

+  i[cos ( 6 )cosh( 6 )-sin( 6 )sinh( 6 ) ]  }.  (B-19) 

This  completes  the  evaluation  of  bracket  number  4 . 

Having  evaluated  each  of  the  inner  brackets  1,  2,  3,  and  4  in 
the  denominator,  equation  B-4,  the  following  relations  have  been 
established: 


i  ),  -  -1  )3  -  X  . 

(B-20) 

(  h  ■  (  )4  -  v  • 

(B-21) 

These 

these 

equations  shall  define  the  two  parameters  X  and 
parameters,  the  denominator  takes  the  following 

y.  Using 
form 

v  =  yHQ(kb)  [yJ0(ka)\  +  kJ£(ka)y]  + 

kH^ (kb)  [kJ^(ka)X  -  yJQ(ka)y] 

(B-22) 

*  X[y*HQ(kb)J0(ka)  +  kaHx (kb) jx (ka ) ]  + 

-yky[H0(kb)J1 (ka)  -  Hx (kb) JQ(ka) ] 

or 

V  =  ka  +  yp 

(B-23 ) 
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where 


a  =  y2H0<kb)J0<ka)  +  k,H1(kb)J1(ka) ,  (B-24) 

p  =  yk[J0(ka)H1(kb)  -  HQ (kb) (ka ) ] .  (B-25) 


Equation  B-23  is  the  final  result  for  the  denominator.  The 
parameters  a  and  p,  given  by  equations  B-24  and  B-25,  contain 
products  of  Bessel  functions  of  real  arguments.  In  general,  these 
must  be  calculated  by  the  methods  described  in  Appendix  A.  The 
denominator  is  calculated  by  first  calculating  X,  y,  a,  and  p, 
from  equations  B-10,  B-13,  B-24,  and  B-25,  respectively,  and 
then  using  B-23. 

The  next  element  in  the  calculation  of  d0  is  the  numerator, 
which  is  given  by  equation  40  with  n=0.  Evaluating  the 
determinant  gives  U«  = 


yk [ HQ ( kb ) Jq ( kb )  -  J0(kb)H^(kb)][J0(ya)H^(ya)  -  HQ (ya) ( ya) ] 


-  yk[H0(kb)J1(kb)  -  J^kbJH^kbjjfj^yajH^ya)  -  H()(ya)J1(ya)  ]; 


which,  using  the  definition  A-5,  becomes 


«  yk[J1(ku)Y0(kb)  -  J0(kb)Y1(kb)][J0(ya)Y1(ya)  -  YQ (ya) J1 (ya) ] . 


The  expressions  in  brackets  may  be  simplified  by  using  the 
identity 


[J 


n+l<z>Vz) 


-  Jn<z>Yn+l 


(z)  ] 


2/ttz  . 


Therefore,  the  numerator  takes  the  final  form 


”4 


4 

7isab  ' 


( B-26 ) 
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This  result  is  exact.  The  coefficient  d0  is  the  quotient  of  the 
numerator  B-26  and  the  denominator  B-23.  In  practice,  however,  it 
is  easier  to  calculate  the  reciprocal  of  d0,  which  is  given  by 


1 

3e 


(Xa  +  jap ) 


7Taab 

4 


(B-27 ) 


It  is  important  to  specify  the  conditions  under  which  this 
result  is  valid.  Essentially,  equation  B-27  is  a  good 
approximation  as  long  as  the  principal  asymptotic  forms  are 
valid.  Generally,  the  error  in  the  principal  asymptotic  forms  for 
the  Bessel  functions  decreases  for  increasing  |z|.  For  example, 

when  | z | > 1 25 ,  the  error  is  less  than  or  equal  to  3xio~J; 

when  | z | >1250 ,  the  error  is  less  than  or  equal  to  3xi0-4;  and 

when  |z | >12500,  the  error  is  less  than  or  equal  to  3xio-5.  It 
will  be  assumed  that  the  results  of  this  appendix  are  valid  when 
| z | >125 .  This  choice  for  the  validity  criterion,  while  somewhat 
arbitrary,  is  substantiated  by  accurate  calculations  of  l/d0 
which  show  that  for  all  the  cases  studied  in  this  report  the 

error  in  B-27  is  less  than  or  equal  to  10-5  when  |z|-l25.  For  the 
arguments  z-ya  and  z=yb,  this  condition  is  satisfied  when 
|ya|>125,  which  implies 


ya  |  =  /coy  eoa  >  1  25 . 


( B— 2  8 ) 


For  the  ranges  of  physical  parameters  considered  in  this  appendix 
this  criterion  places  a  constraint  on  the  minimum  frequency  f 
for  which  equation  B-27  is  applicable.  For  a  given  conductivity 
o,  and  a  given  radius  r=a.  Table  B-l  gives  the  minimum  frequency 
obtained  by  equation  B-28. 

Finally,  the  resonance  condition  for  the  TM  modes  will  be 
derived.  It  is  known  that  electromagnetic  resonance  occurs  for 
certain  discrete  frequencies,  and  that  this  causes  a  reduction  in 
the  shielding  effectiveness  at  those  frequencies.  The  shielding 
effectiveness  may  be  reduced  by  3  or  4  orders  of  magnitude  due  to 
this  effect.  By  comparison  with  the  resonance  condition  for  a 
pair  of  parallel  conducting  planes  separated  by  a  distance  2a,  it 
is  logical  to  expect  a  resonance  to  occur  when  the  inner  diameter 
of  the  cylinder  is  an  integral  number  of  half 
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Table  B-l .  Minimum  frequency  for  a  given  cylinder  radius. 


a  mm 

—.myji 

f  Hz 

63.0 

lxlO7 

£  5xl04 

44.4 

lxlO7 

£  1x10* 

19.9 

lxlO7 

^  5x10 8 

14.0 

lxlO7 

^  1x10* 

4.45 

lxlO7 

_ 

£  lxlO7 

wavelengths,  that  is,  2a  =  n(X/2).  Using  the  result  for  the 
shielding  effectiveness ,  equation  B-27,  it  may  be  shown  that  the 
resonance  condition  for  a  pair  of  concentric  cylinders  is 
J0(ka)  *  0.  To  prove  this,  first  note  that  as  the  frequency 
increases  from  zero,  the  first  resonance  occurs  approximately 
when  (X/2 )~2a,  or  f~(c/4a).  For  such  large  frequencies,  that  is, 
f > (c/4a) , 


6  =  (“|ils.)1/2(b-a)  >  (TT-C°^a)I/2A  >  4.7  ,  ( B-29 ) 


where  it  is  assumed  that  A  =  (b-a)/a  >  0.001,  a>7.5  mm,  and 

o> 10 7  ohm-1m-1.  This  frequency  range  will  be  called  the  resonance 

regime.  If  terms  of  the  form  e-6  are  neglected  in  the  defining 
relations  B-20  and  B-21,  for  X  and  y,  respectively,  then  this 
shows  that  y=iX.  Therefore,  the  inverse  shielding  effectiveness 
is  approximately  proportional  to 


|Xa  +  y(3 1  -  |X(a+iP)  j  =  |X||a+ip|.  (B-30) 

Since  |X|  -  e6  is  an  increasing  function  of  frequency,  the 
resonance  phenomena  must  arise  from  the  a+ip  term.  Studying  this 
term  in  detail  shows  that  since,  for  frequencies  well  below  the 
infrared,  |y’|'  |yk  |  >  >  |  k*  |  ,  ':he  dominant  term  in  a+ip  is 
the  y*  term: 
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a+ip 


y*H0 (kb) J0 (ka) . 


(B-31 ) 


Consequently,  the  minima  in  a+ip  occur  when  J0(ka)  «  0.  This 
establishes  the  resonance  condition  for  the  TM  modes.  Note  that 
since  the  zeros  of  J0(x)  and  J,(x)  alternate  for  x  real,  neither 
a  or  p  is  ever  zero,  and  consequently  the  shielding  effectiveness 
never  goes  to  zero.  The  resonant  frequencies  are  given  by 


ka 

*  j  o (n)  , 

n  *  1,2,3, .. .  ; 

(B-32) 

or 

fn 

-  2iS  ' 

n  =  1,2,3,...  ; 

( B—  33) 

where  j0(n)  are  the  nth  zeros  of  J0(x). 
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APPENDIX  C 
CALCULATION  OF  D0 
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The  calculation  of  D0  is  similar  to  the  calculation  of  d0 
given  in  Appendix  B.  The  magnitude  of  D0  is  essentially  the 
inverse  shielding  effectiveness  in  the  TE  case,  as  |d6 |  was  in 
the  TM  case.  D0  is  given  by  equation  54.  By  expanding  the 
determinant  in  equation  59,  the  denominator  is  given  by  Y= 


H0<kb>l^l5Fl  {  [Ji<Yb>H'(ya)  -  JJ(va)H’  (yb)  ]j  + 

-  ll?„j6(ka)  [J0(Ya>H6<Yb)  -  Ji(Yb)H0(ya)]2)  + 

-llE„Ho(kb>  (  rik„J6(ka>  tJ0,ya)H0(Yb)  -  J0,Yb)H0(ya)]3  + 

J0(ka>loic)  [J0(Yb)H^(ya)  -  J^(ya)H0(yb)  ]4}  .  (C-l) 

The  inner  brackets  are  numbered  1,2,3,  and  4  for  easy  reference. 
Note  that  these  brackets  are  equivalent  to  the  brackets  1,2,3, 
and  4  in  Appendix  B.  Using  the  notation  from  Appendix  B,  the 
expression  for  the  denominator  C-l  may  be  written 

H0 (kb)  U-icoe )  (J0 (ka)  ^o-iwe )  l  )l  "  IwE0J6(ka)  (  ^  2^  + 

-ife.H6<kb’  tife„Ji<ka>  (  )3  -  Jo<ka>(dfei  !  >4  )•  <c-2> 


Using  the  principal  asymptotic  forms  B-2  and  B-3  to  evaluate  the 
inner  brackets,  exactly  as  was  done  in  Appendix  B,  leads  to 


( )l 

-  -( 

)3  *  X  , 

(C-3 ) 

( )2 

-  ( 

)4  -  v  ; 

( C-4 ) 

which  are  identical  to 
calculations  are  given 

B-20  and  B-21 .  The  details 
ii.  Appendix  B.  Using  these 

of  these 
results  the 
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denominator  becomes  Y= 


<orfc)  Vkb> 

{ 

f  — 1—  ) 
lO-iWE ' 

XJQ(ka) 

k 

icoE  o 

MJ^fka) 

}  ♦ 

ISH.  Ho(kb| 

{ 

f  V,  1 

lO-iW£ 1 

pJ0(ka) 

+ 

iWE0 

Xj^(ka) 

)  . 

(C-5) 

which  may  be  written 


Y  «  X£  +  yn  / 


(C-6 ) 


where 

5  =  [psHQ(kb)J0(ka)  +  q*H£(kb)J£(ka)] 

-  [pIH0(kb)J0(ka)  +  q*H1(kb)J1(ka)]  ,  (C-7) 

H  -  pq  [J0(ka)H^(kb)  -  Jq (ka)HQ (kb) ] 

-  pq  [j1(ka)H£)(kb)  -  JQ (ka)H1  (kb)  ]  ,  (C-8) 

and ,  where 

p  -  -  ana  9  -  is.-  'c-9> 

This  completes  the  evaluation  of  the  denominator.  The  evaluation 
of  the  numerator,  which  is  given  by  the  determinant  in  equation 
58,  procedes  exactly  as  in  Appendix  B.  The  numerator  may  be 
calculated  exactly  with  the  result 


Tikb*  *  ^  Trya*-i)  *  “yf  n*ab  =  ia>E0(o-i«£)  nJab  10^ 

The  coefficient  D0  is  the  quotient  of  the  numerator  and  the 
denominator  given  by  equations  C-6  and  C  10,  respectively.  In 
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practice  it  is  more  convenient  to  work  with  the  inverse  of  D„ . 
The  inverse  of  D0  takes  the  form 


1 

Do 


(Xa'+  yp' ) 


TT*ab 

4 


where 

Q'  "  "pq  C  *  -yk[(p/q)H0(kb)J0(ka) 

P'  -  n  =  -ykfj^kajH^kb) 


-  (q/pJH^kbjJ^ka)  ] 

-  J0(ka)H1(kb)]. 


(C-ll) 

( C-12 ) 
(C-13) 


A  further  simplification  may  be  made  by  noting  that 


(q/p) 


r o-icoe  ^  k 
'  i«E0  ‘  y 


Therefore 

a'  -  k2H0(kb)JQ{ka)  +  y *1^  (kb) Jx (ka) 
p'  =  yk[j0(ka)H1(kb)  -  J2 (ka)HQ(kb) ] 


( C— 14) 

(C-15) 
(C-16 ) 


The  result  for  D0 ,  equation  C-ll,  is  identical  to  the  result  for 
d0,  equation  B-27,  except  that  a'  and  P'  in  the  former  equation 
replace  a  and  p  in  the  latter  equation. 

As  in  Appendix  B,  it  is  important  to  specify  the  conditions 
under  which  this  result  is  valid.  The  equation  for  De,  equation 
C-ll,  will  be  assumed  valid  for  | ya  j  >125.  Just  as  in  Appendix  B, 
for  a  given  radius  r-a  and  a  given  conductivity  o,  this  imposes  a 
restriction  on  the  minimum  frequency  for  which  C-ll  is  valid. 

This  leads  again  to  the  inequality  /cooy0a>125.  and  the 
corresponding  Table,  Table  B-l. 

Lastly,  the  resonance  condition  for  the  TE  modes  will  be 
derived.  As  in  the  TM  case,  electromagnetic  resonance  will  occur 
for  certain  discrete  frequencies,  causing  a  reduction  in  the 
shielding  effectiveness  at  tho^e  frequencies.  As  in  the  TM  case, 
the  shielding  effectiveness  may  be  reduced  by  3  or  4  orders  of 
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magnitude  due  to  the  resonance  effect.  Using  the  result  for  the 
shielding  effectiveness,  equation  C-ll,  it  may  be  shown  that  the 
resonance  condition  is  J^ka)  ■  0.  The  derivation  of  this  result 
is  exactly  the  same  as  for  the  TM  case  in  Appendix  B.  Based  on 
the  resonance  condition  for  a  pair  of  parallel  conducting  planes 
it  is  known  that  the  first  resonance  occurs  approximately  when 
(X/2)~2a,  or  f~(c/4a).  For  large  frequencies  such  that  f>(c/4a), 
it  may  be  assumed  that  6 >4. 7,  as  was  done  in  Appendix  B,  equation 

B-29 .  Therefore,  in  the  resonance  regime,  terms  of  the  form  e-6 
may  be  neglected  in  the  expressions  for  both  X  and  y.  This  leads 
to  the  result  y=iX.  Therefore  the  inverse  shielding  effectiveness 
is  proportional  to 


Xa'  +  yp ' |  =  |X(a'+ip')|  «  jxjja'+ip'j.  (C-17) 


r 

Noting  that  |X|  -  e  is  an  increasing  function  of  frequemcy,  the 
resonance  phenomena  must  arise  from  the  a'+ip'  term.  Since,  for 
frequencies  well  below  the  infrared,  | y 2 [ > > | yk | > > | k * | ,  the 
dominant  term  in  a'+ip'  is  the  y*  term: 


a'+ip*  -  y*^  (kb)J2  (ka) .  (C-18) 


Consequently,  the  minima  in  a'+ip'  occur  when  J,(ka)  -  0.  This 
establishes  the  resonant  condition  for  the  TE  modes.  Note  that 
since  the  zeros  of  J0(x)  and  J,(x)  alternate  for  x  real,  neither 
a'  or  p'  is  ever  zero,  and  consequently  the  shielding 
effectiveness  is  never  zero.  The  resonant  frequencies  are  given 
by 


ka 

-  j.(n)  , 

n  *  1,2,3,...  ; 

(C-19 ) 

or 

fn 

“  2tt£  ^  1  (n)  ' 

n  -  1,2,3,...  ; 

(C-20 ) 

where  j»(n)  are  the  nth  zeros  of  J,(x). 
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APPENDIX  D 

APPROXIMATE  FORMULAS  FOR  SHIELDING  EFFECTIVENESS 


66 


It  is  extremely  useful  to  have  an  approximate  closed  form 
result  for  the  shielding  effectiveness.  Such  a  form  is  easily 
derived  in  the  TM  case  by  using  the  results  of  Appendix  B.  The 
shielding  effectiveness,  given  by  equation  B-27,  is 


SE~1/2  *  j\a  +  yp|,  (D-l ) 


which,  for  frequencies  in  the  resonance  regime,  may  be  written 


SE  1/2  «  | X |  |a+ip  |  .  ( D-2 ) 


Now  substitute  the  result  from  Appendix  B, 


M  =  i  M'1  (ab)"1/2  e6,  (D-3 ) 

which  is  valid  in  the  resonance  regime;  and  based  on  the  result 
B-31 ,  make  the  approximation 


|a+ip|  -  |y2  |  M0  (kb)Mq  (ka )  -  | y*  |  (^) 1/2 (^) 1/2 ,  (D-4) 


where,  for  x  real, 

M0(x)  -  | H0 ( x ) |  *  (J2(x)  +  Y2(x)j1/2,  (D-5) 


and  where  it  is  assumed  that 


I J0 ( X ) I  ~  M0(x). 

This  yields  the  simple  formula 


SE'172 


1_  t  a  >  1/2 
2  WE  o 


(D-6) 


(D-7) 
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or. 


where 


SE 


4 COE  a  e-26 


1/2 


(b-a) . 


(D-8) 


(D-9) 


This  shows  that  the  shielding  effectiveness  is  independent  of  the 
radius  of  the  cylinder  and  depends  only  on  the  thickness  of  the 
shield.  A  formula  for  the  shielding  effectiveness  in  dB  is: 


SEdB  "  io{iog( 


4E , 


•)  -  log(w)  +  (2oy0 )1/2(b-a)log(e)  «1/2}.  (D-10) 


This  formula  is  a  good  approximation  in  the  resonance  regime, 
that  is,  f>(c/4a).  However,  the  approximation  breaks  down  for  low 
frequencies.  If  equation  D-8  is  considered  as  a  function  of  to, 
that  is,  f(w) ,  then  this  function  has  a  minimum  at  the  frequency 


«0 


2  1 
y0o  (b-a) 2 • 


(D-ll) 


As  «  goes  to  zero,  f(«)  approaches  infinity.  Therefore,  the 
formulas  D-8  or  D-lo  should  not  be  applied  for  frequencies  less 
than  w0.  The  error  at  this  frequency  is  roughly  38%  for  the  cases 
studied  in  this  report.  Consequently,  this  formula  only  provides 
a  rough  estimate  in  this  lower  frequency  range. 

The  same  formulas  for  the  shielding  effectiveness  may  be 
derived  for  the  TE  case  by  using  the  results  of  Appendix  C.  The 
shielding  effectiveness,  given  by  equation  C-ll,  is 


which. 


SE_1/2 
for  frequencies 


«  Hlab  |Xa,+  m(3,  |  #  (D-12 ) 

in  the  resonance  regime,  may  be  written 
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|x|  ja'+ip'  | 


( D-l 3 ) 


SE~1/2  = 


Now  make  the  approximations  given  by  D-3,  D-4,  D-5,  and  D-6 
again  yields  the  formulas  given  in  equations  D-8  and  D-10. 


.  This 
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APPENDIX  E 

AN  ALTERNATE  DEFINITION  OF  SHIELDING  EFFECTIVENESS 
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An  alternate  definition  of  shielding  effectiveness  may  be 
given  in  terms  of  the  energy  density  of  the  incident  and  axial 
fields.  The  most  natural  definition  of  shielding  effectiveness  is 


SE 


(2) 


U 


axis 


u 


me 


(E-l) 


where  U  is  the  energy  density  of  the  fields.  The  factor  of  two  is 
included  so  that  the  shielding  effectiveness  approaches  1  as  the 
frequency  approaches  zero.  The  average  energy  density  of  an 
electromagnetic  field  is  defined  by  the  relation 

U  =  •|{E»D  +  H»B},  (E-2 ) 

where  E,  D,  H,  and  B  are  the  real  fields,  and  {}  denotes  a  time 
average.  In  terms  of  the  complex  fields,  and  in  free  space,  this 
is  equivalent  to 


U  -  i[e0E*E*+  y0H*H*].  (E-3 ) 


It  is  interesting  to  compare  the  definition  E-l  to  the 
definition  used  in  Section  7  for  the  TM  and  TE  modes,  equations 
140  and  141.  It  will  be  shown  that  these  definitions  are 
equivalent.  The  energy  density  of  the  incident  planewave  is 


U 


inc 


( E-4  ) 


For  the  TM  case  the  E-field  at  r=0  is  given  by  the  n=0  term  in 
the  series  expansion  since  Jn(0)-0  for  n=±l,±2,...  .  The  4>- 

component  of  the  H-field  is  given  by  the  n=±l  terms  since 

otherwise.  The  energy  density  of  the  field  on  the  axis  of 

the  cylinder  is 

(TM)  Uaxis  =  ^e0e2  [|dQ|2+  |d_1+  d+1 1 2  ]  .  (E-5 ) 
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Similarly,  for  the  TE  case  the  H-field  at  r-0  is  given  by  the 
n-0  term  in  the  series,  and  the  ^-component  of  the  E-field  is 
given  by  the  n  *  ±1  terms.  Therefore  the  energy  density  on  the 
cylinder  axis  is 


(TE) 


U 


axis 


UDol 


d-i+ 


+1 


I2] 


(E-6 ) 


The  symmetry  about  the  x-z  plane  implies  that  the  fields  are 
represented  by  a  cosine  series.  This  implies  that  cl_1«  -d^  and 

D_1-  -D^  Therefore,  the  shielding  effectiveness  for  the  TM  and 

TE  modes  are  given  by 


(TM) 

SE 

-  Kl2 

(E-7  ) 

(TE) 

SE 

-  lDol2- 

(E-8 ) 

In  terms  of  dB: 

(TM) 

SEdB  - 

-20  log 

ldo  1  ' 

(E-9) 

and 

(TE) 

SEdB  = 

-20  log 

l°0l  • 

( E-10 ) 

These  formula  are  equivalent  to  equations  146  and  147.  Therefore 
the  definitions  140  and  141  given  in  section  7  for  the  shielding 
effectiveness  of  the  TM  and  TE  modes,  respectively,  are 
equivalent  to  E-l. 

Next,  consider  a  perpendicularly  incident  wave  with  an 
arbitrary  state  of  elliptical  polarization.  The  definition  E-l 
lends  itself  naturally  to  such  a  situation.  It  will  be  shown  that 
the  approximate  formula  D-7  for  the  shielding  effectiveness  is 
true  for  any  state  of  polarization  of  the  incident  wave.  To 
begin,  note  that  the  TM  incident  wave  has  rectangular  components 
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E2  =  E« 


(E-ll) 


Hy  -  -Ho, 

and  the  TE  incident  wave  has  the  components 


(E-12) 


H2  "  H» 


( E— 13) 


Ey  -  E0. 


(E-14 ) 


where  E„/H0=  J\}0/e0  .  The  electromagnetic  fields  for  an  incident 
planewave  having  an  arbitrary  state  of  elliptical  polarization 
may  be  written 


c1{tm)0+  c2{te}0 


(E-15 ) 


where  the  bracket  {}Q  stands  for  the  totality  of  field 

components,  that  is,  both  E  and  H,  for  either  the  TM  or  TE 
incident  waves.  The  coefficients  c,  and  c,  are  complex  constants 
such  that 


|c,|2+  |c,|2  >  1. 


(E-16) 


Using  the  explicit  TM  and  TE  field  components  given  above,  the 
incident  wave  may  be  written  in  the  form  ! 


Einc  '  E"(c‘ez+  <=.«v) 


( E— 17) 


Ho(c,e  +  c,e  ) . 


2  1  y 


(E-18) 


where  e„  and  e„  are  the  unit  vectors  in  the  y  and  z  directions, 
y  z  J 

respectively.  The  energy  density  of  the  incident  wave  ,1s  computed 

using  equation  E-3,  which  yields 
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ulnc  ■  ^.|Et|2(|c,|2+  |c,|2)  -  jE(  I  E„  | 2  •  (E-19) 

Now  energy  density  of  the  axial  fields  must  be  calculated.  If  the 
incident  wave  is  of  the  form  E-15  then  by  the  principle  of 
superposition  the  solutions  for  the  total  fields  have  the  form 


c, {tm}  +  c,{te) , 


(E-20) 


where  inside  the  cylinder,  for  example,  {tm)  stands  for  the 
totality  of  field  components,  both  E  and  H,  given  by  equations 
29,  30,  and  31;  and  {te}  stands  for  the  totality  of  field 
components  given  by  equations  47,  48,  and  49.  Using  this 
representation  the  fields  on  axis  are  given  by 


Ez  =  CjE0d0,  (E-21 ) 

E  =  c2E0D0 .  (E-22) 

Ct 


Consequently,  the  energy  density  on  the  axis  of  the  cylinder  is 


“axis*  1  |E0 | 2(tCl | 2 |d„ | 2+  |c,|2|d„|2}. 


(E-23 ) 


Therefore,  the  shielding  effectiveness  for  a  perpendicularly 
incident  planewave  having  an  arbitrary  state  o.f  elliptical 
polarization  is 


SE  - 


(E-24) 


As  expected,  this  is  a  superposition  of  the  shielding 
effectiveness  for  the  independent  TM  and  TE  modes.  When  the 
approximate  formula  D-8  of  Appendix  D  is  valid,  then  |d0|-|D0| 
which  implies  that  E-24  takes  the  form 


SE  =  |d0 |2  -  | D0 | 2 .  ( E-25 ) 
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This  proves  that  the  approximate  formula  D-8  is  valid  for  any 
state  of  polarization  of  the  incident  wave. 
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APPENDIX  F 

APPROXIMATIONS  FOR  THE  SPHERICAL  BESSEL  FUNCTIONS 
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As  was  shown  in  section  7,  the  fields  at  the  origin  are 
determined  by  the  two  coefficients  p,  and  q, .  Therefore  the 
calculation  of  the  shielding  effectiveness  is  concerned  primarily 
with  the  numerical  evaluation  of  these  two  quantities.  According 
to  the  equations  for  p,  and  qw  equations  209  and  212,  the  main 
problem  is  to  evaluate  the  determinants  W,  and  W,  given  by 
equations  205  and  206.  To  do  this,  special  care  must  be  taken 
when  computing  each  of  the  four  brackets  [ ] 1 ,  []2,  []3,  and  []4« 

By  using  the  Rayleigh  formulas,  the  spherical  Bessel  and 
Hankel  functions  and  the  required  derivatives  are  given  by 


j1(z)  = 

sinz  cosz 

za  z  ' 

(F-l) 

[zj1(z>]'  - 

(l  -  i,)sinz  +  ^-cosz. 

( F-2 ) 

h1 ( z )  = 

-Mi.-  |)eiz. 

(F-3) 

[zh1(z)]‘  = 

M|.-  |  -  Deiz. 

(F-4) 

In  general,  these  formulas  may  be  used  to  calculate  any  one  of 
these  quantities  exactly.  However,  for  the  complex  arguments 
z«ya  and  z~yb,  it  is  not  practical  to  compute  these  quantities 
directly  from  the  above  equations,  except  at  low  frequencies.  The 
reason  is  that  for  high  frequencies  the  values  of  the  complex 
arguments  z*ya  and  z=yb  yield  values  of  the  spherical  Bessel 
function  j  A  <  z )  which  are  astronomically  large  since  they  have  a* 

large  positive  exponent,  and  values  of  the  Hankel  function  h^z) 

which  are  astronomically  small  since  they  have  a  large  negative 
exponent.  The  magnitudes  of  these  numbers  is  much  greater  than 

the  range  of  most  computers;  on  the  order  of  io10'®0®.  However, 
the  Bessel  and  Hankel  functions  with  complex  arguments  have  been 
grouped  together  in  the  brackets  []1#  []2,  []3,  and  []4,  in 

equations  205  and  206,  and  since  products  of  these  functions 
appear  in  pairs  the  large  positive  and  negative  exponents  cancel. 
By  making  a  simple  approximation  the  products  of  the  Bessel  and 
Hankel  functions  may  be  calculated  quite  easily  and  accurately. 
First  note  that  for  a  good  conductor,  and  for  frequencies  such 
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that  (o/«e0)<<1,  the  propagation  constant  inside  the  conducting 
medium,  y,  is  given  by 


V  -  -  filial  -  I!  UjJt). 


(F-5) 


which  implies  the  approximation 


v  -  5  (ff)l/2  ^ 


(«MeO) 


il 

1/2  14 

e  , 


(F-6 ) 


Considering  the  following  ranges  of  physical  parameters: 

a  >  7.5  mm,  f  >  5*104  Hz,  and  o  >  10*  ohm-1m-1,  implies  that  the 
arguments  of  the  spherical  Bessel  functions  satisfy  the 
following  inequalities:  ka  >  0.139  and  |z|  -  |ya|  >  14.89. 
Therefore,  writing  z  in  the  form  z  *  £+i£,  these  conditions 

imply  £  >  10.0.  Consequently,  if  terms  of  the  form  e^Z“  e^e~^ 

-iz  iE  +E 

are  neglected  compared  to  the  terms  of  the  form  e  *  e 

then  the  formulas  F-l  and  F-2  may  be  approximated  by  the  formulas 


j^z)  *  §(§,+  j)e~iZ' 

[ z j ! ( z ) ] 1  -  -f(|a  +  I  -  1 ]e_iZ • 


(F-7) 


(F-8  ) 


These  approximations  shall  be  used  to  evaluate  the  brackets  in 
equations  205  and  206.  Proceeding  with  the  evaluation  of  bracket 
number  1,  let  z,*ya  and  z2=yb,  then  using  F-2,  F-4,  F-7,  and  F-8 
yields 


Fj(zl)r^(za)  —  rj (za )Th(z, ) 


(F-9 ) 


[l  -  z72-  zl2+  (zlzJ)"A+  (Z, Za)“z]  isin(z*-Zj )  + 


-  [zT1-  z7x+  z72z71-  z71z72]  icos(z,-z,). 


(F-10 ) 
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Similarly,  simplifying  bracket  number  two  yields 


[J2  -  h 1(2l)rj(z4)  -  j1(zI)rh(za)  ( F-ll ) 

*  -[zi  ~  (Z[ z2 )  —  (z,z2)  ]  isin(za— z, ) 

-  [z71+  z72z;1-  z71z72]  icos(za-z,).  (F-12 ) 

Bracket  number  three  is  given  by 

[ ] 3  -  rh(zl)j1(za)  -  rj(z1)h1(zI)  ( F— 13) 

*  _[z j  ~  (ZiZj )  —  (Z|Z2)  ]  isin(za-z , )  + 

+  [zl1-  z72z71+  z71z72]  icos(za-z,).  ( F-14 ) 

And,  bracket  number  4  is  given  by 

[]4  =  j1(zl)h1(za)  -  h1(z»)j1(za)  (F-15 ) 

-  [ ( z i z a ) ~  +  (Z)Zj ) -  ]  isin(za-Zj )  + 

\ 

—  0  — 1  — 1  —  o 

-  [z,  Za  -  z7  z2  ]  icos(za-z»).  ( F-16 ) 


In  general,  these  approximations  are  very  good  for  |ya| >14.14. 
These  formulas  for  each  of  the  four  brackets  are  used  in  the 
calculation  of  the  shielding  effectiveness  of  a  spherical 
conducting  shell  at  high  frequencies.  It  is  not  necessary  to 
separate  the  real  and  imaginary  parts  of  these  equations  since 
all  numerical  work  is  performed  on  a  computer  which  supports 
complex  arithmetic. 

In  addition  to  the  approximations  F-7  and  F-8,  for  the 
spherical  Bessel  functions,  another  set  of  approximations  is 
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required  when  the  magnitudes  of  the  arguments  is  very  small, 
that  is,  |z | <<1.  Approximations  of  this  type  are  needed  to 
calculate  the  shielding  effectiveness  for  very  low  frequencies 
where  ka<<l.  The  problem  at  very  low  frequencies  is  not  with  the 
evaluation  of  the  bracketed  terms,  but  with  the  coefficients  of 
these  brackets  in  equations  205  and  206.  These  coefficients  are 
products  of  spherical  Bessel  functions  having  the  real  arguments 
x*ka  and  x*kb.  If  x<<1,  problems  arise  if  j^x)  and  [xj^x)]’  are 

computed  by  using  equations  F-l  and  F-2  on  a  calculator  or 
computer  which  has  an  accuracy  of  12  significant  figures.  For 

example,  if  x«10-6,  then  on  a  machine  with  an  accuracy  of  12 
significant  figures,  a  simple  algorithm  based  on  F-l  and  F-2 
yields  j1(x)-=0  and  [xj  1  (x)  J '  =0 .  This  situation  is  not  acceptable 

since  it  may  lead  to  erroneous  results  for  the  shielding 
effectiveness.  Therefore,  for  |z|<<1,  the  Bessel  and  Hankel 
functions  shall  be  approximated  by  the  relations 


j  i  (2)  * 

( F-l 7 ) 

[zj  » (z) ]'  * 

¥■ 

(F-18 ) 

h,  (z)  = 

3*  -  M|.+  f]. 

(F-19 ) 

[zh» (2) ] '  * 

¥ +  Mp-  ¥■ 

( F-20 ) 

These  approximations  may  be  obtained  from  the  defining  relations 
F-l,  F-2,  F-3 ,  and  F-4  by  straightforward  series  expansions  of 
the  trigonometric  and  exponential  functions.  These  equations  are 

utilized  in  numerical  computations  when  |z|<10  . 

The  resonance  conditions  shall  now  be  derived  for  the  electric 
and  magnetic  fields.  By  analogy  with  the  interference  effects  in 
a  Fabry-Perot  interferometer  it  is  expected  that  the  fields 
inside  the  spherical  cavity  should  exibit  resonances  for  certain 
discrete  frequencies.  Consider  the  resonances  in  a  system  of  two 
parallel  planes  seperated  by  a  distance  d.  Resonances  will  occur 
when  the  separation  of  the  planes  is  equal  to  an  integral  number 
of  half-wavelengths,  that  is,  n(X/2)«d,  or  f=(c/2d)n.  Therefore, 
it  is  logical  to  expect  that  the  first  resonance  for  the  sphere 
will  occur  near  the  frequency  f=(c/4a),  or  equivalently 
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ka= (n/2 )*1 .  For  frequencies  such  that  f£(c/4a),  the  magnitude  of 
the  arguments  z=ya  satisfies  the  inequality 


|  ya  |  =  (woy0]1/2a  *  (Ttacoy^l/2.  (F-21) 


7  -1-1 

For  the  parameters  o=10  ohm  m  and  a=7.5  mm,  this  inequality 
yields  | ya | >6664.  Therefore,  in  the  resonance  regime,  the 
approximations  F-10,  F-12,  F-14,  F-16  will  certainly  hold  for  the 
brackets  []1#  [)2,  []3,  and  []4. 

The  resonances  for  the  magnetic  field  shall  occur  when  the 
magnitude  of  the  coefficient  p,  is  a  maximum,  since  at  the  center 
of  the  sphere  | H/H0 1  =  I p » | .  This  implies  that  the  reciprocal  of 
|p, |  will  be  a  minimum.  To  first  order. 


W,  =  j^kajh^kb)  [  ]1  -  j1(ka)h1(kb)  isin[y(b-a)  ] ,  (F-22) 


according  to  F-10,  F-12,  F-14,  and  F-16.  Therefore  the  resonances 
for  the  magnetic  field  will  occur  when  the  quantity 

I P 1  I ~ 1  -  | ya | (kb) | w, |  -  |ya| (kb) |31(ka) I |h1(kb) I |sin[y(b-a)] I 

(F-23) 

is  a  minimum.  Clearly,  this  occurs  when  j1(ka)=0.  This  is  the 

resonance  condition  for  the  magnetic  field.  Note  that  this  is 

also  the  resonance  condition  for  the  magnetic  shielding 

effectiveness  SE  .  since  SE  =  | p .  1 3 . 

m  m  1  c  1 

The  resonances  for  the  electric  field  shall  occur  when  the 
magnitude  of  the  coefficient  q,  is  a  maximum.  This  implies  that 
the  reciprocal  of  |q, |  will  be  a  minimum.  To  first  order. 


W a  =  l\(ka)rh(kb)jj£  []4 

-  (ka)rh(kb){^-[  (ya)  (yb)  ]_1isin[y (b-a )  ] ,  (F-24) 
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according  to  F-10,  F-12,  F-14,  and  F-16.  Therefore  the  resonances 
for  the  electric  field  will  occur  when  the  quantity 


|q, )  1  =  | k/y | (ka) (kb) |w, | 

-  |y/k||rj(ka)||rh(kb)||sin[y(b-a)]|  (F-25) 


is  a  minimum.  Clearly,  this  occurs  when  r\j(ka)=0,  that  is, 

[ p J ^  (p)  ]p.jca“0 .  This  is  the  resonance  condition  for  the  electric 

field.  This  is  also  the  resonance  condition  for  the  electric 
shielding  effectiveness  SE  ,  since  SE  =  | q ! | * . 

For  the  electromagnetic  shielding  effectiveness,  which  is 
defined  in  terms  of  the  energy  density  of  the  fields,  the 
resonances  will  occur  when  either  of  the  previously  derived 
conditions  are  satified  separately  ,  since  SE  =  SEm+  SEe. 

Therefore,  the  resonance  conditions  are 


j^ka)  - 
[ Pj 1 ( P) ]p=ka~ 


0,  (F-26) 

0.  ( F-27 ) 


The  first  few  zeros  of  F-27  are  given  in  Table  F-l. 


Lastly,  an  approximate  formula  for  the  SE  of  a  sphere  shall  be 
derived  which  is  identical  to  the  approximate  formula  for  the 
cylinder.  This  approximate  formula  shall  be  derived  for 
frequencies  in  the  resonance  regime.  For  the  magnetic  shielding 
effectiveness,  start  with  equation  F-23: 

|p. r1  *  | ya | (kb) |Wj |  -  |ya| (kb) I jx(ka) I |h1(kb) I |sin[y(b-a) ] I 

(F-28 ) 

Assuming  that  |y (b-a) | >14 .14,  the  sine  term  may  be  approximated 
by  a  single  exponential. 
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n 

zn 

1 

2.743707  269992 

2 

6.116764  264462 

3 

9.316615  628566 

4 

12.485937  368200 

5 

15.643866  106348 

6 

18.796253  353454 

Table  F-l .  Zeros  of  [zj1(z)]'  for  z  real,  z>0. 


sin[ y ( b-a ) ] j 


1 

2 


Y  |  2 


-1/2 


(b-a) 


(F-29) 


The  magnitude  of  the  spherical  Bessel  and  Hankel  functions  may  be 
roughly  approximated  by  their  modulus  ([19],  page  439,  equations 
10.1.28  and  10.1.29) 


ji(x) 

hx(x) 


Assuming  ka>>l,  then  to  first  order  j 1 (ka ) = ( 1/ka )  and 
h1(ka)=(l/ka) ,  which  imply  that  equation  F-28  takes  the  form 


-1 


I  Y/k  I 


-1/2 


(b-a) 


( F-30 ) 


This  yields  the  desired  formula  for  the  shielding  effectiveness 
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SE 


m 


4coe  q  e-/2cooy0  (b-a) 


(F-31) 


Next,  the  same  formula  shall  be  derived  for  the  electric 
shielding  effectiveness.  Start  with  equation  F-25: 


|q.|  1  -  | y/k |  | r_j  (ka)  j  | r^(kb)  |  ( sin[ y ( b-a )  ]  |  (F-32) 


Assuming  that  |y(b-a) | >14.14,  the  sine  term  may  again  be 
approximated  by  a  single  exponential  as  in  equation  F-29. 

Consider  the  terms  r\(ka)  and  r^(kb).  Using  a  standard  Bessel 

function  identity,  one  may  write 

[pj1(P)]t  =  pj{(P)  +  j x < P)  «  Pj0(P)  -  jx(P).  (F-33 ) 

[  Ph2  ( p )  ] '  =  phj(p)  +  h2(p)  -  ph0(p)  -  h^p).  ( F-  34 ) 

Then,  roughly  approximating  jQ(p)  and  j1(P)  by  their  modulus,  to 


first  order,  1/p,  since  ka>>l. 

this  shows 

that 

[  P  j  !  ( P )  ] ' 

-  1  -  5  ’ 

1, 

( F—  35) 

[ ph1 <  p ) ] ' 

-  • 

1.  • 

(F-36 ) 

Employing  these  rough  order  of  magnitude  approximations  one 
obtains 


|y/k |  i  «|Y|2'1/2(b-»>, 


which  again  yields  the  formula 


(F-37 ) 
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( F-  38 ) 


SEe 


qj2=  4(oe  a  c-/2<aou  B  ( b-a )  . 


Therefore,  the  final  formula  for  the  electromagnetic  shielding 
effectiveness  in  the  resonance  regime  is 


SE  -  SE  +  SE  -  ^—2-  e“^2“oy°  ( b-a )  . 

me  o 


This  shows  that  the  shielding  effectiveness  at  the  center  of  the 
sphere  is  independent  of  the  radius  of  the  sphere  and  depends 
only  on  the  thickness  and  conductivity  of  the  shield.  In  terms  of 
dB,  this  formula  may  be  written 


SE  (dB)  =  l0{Log(o/8s0 )  -  Log(w)  +  /2 op0 (b-a)Log(e) 


<o1/2}. 


(F-40) 


In  general,  the  formulas  F-39  and  F-40  are  good  approximations 
when  the  conditions  |ya| >14.14  and  | y (b-a ) | >14.14  are  both 
satisfied  simultaneously.  This  implies  that  |ya| >14.14,  or 

1/2 

equivalently  (cooy0)  a  >14.14.  These  formulas  break  down  for  low 
frequencies  and  cannot  be  used  for  frequencies  less  than 


w0 


2  1 
UoO  (b-a) 2  • 
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